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Abstract 

The  paper  presents  a  summary  of  recent  Computational  Fluid  Dynamics  (CFD)  simula¬ 
tions  of  shock  wave  turbulent  boundary  layer  interactions.  This  survey  was  prepared  as  part 
of  the  activity  of  NATO  RTO  Working  Group  10  which  was  established  in  December  1998, 
and  considers  results  obtained  subsequent  to  the  previous  survey  paper  on  the  same  topic  by 
Knight  and  Degrez  (“Shock  Wave  Boundary  Layer  Interactions  in  High  Mach  Number  Flows 
-  A  Critical  Survey  of  Current  CFD  Prediction  Capabilities”,  AGARD  Advisory  Report 
AR-319,  Volume  II,  December  1998).  Five  configurations  are  considered:  2-D  compression 
corner,  2-D  shock  impingement,  2-D  expansion-compression  corner,  3-D  single  fin  and  3-D 
double  fin.  Recent  Direct  Numerical  Simulations  (DNS),  Large  Eddy  Simulations  (LES)  and 
Reynolds-averaged  Navier-Stokes  (RANS)  simulations  are  compared  with  experiment.  The 
capabilities  and  limitations  are  described,  and  future  research  needs  identified. 
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1  Introduction 

Effective  design  of  supersonic  and  hypersonic  air  vehicles  requires  accurate  simulation  meth¬ 
ods  for  predicting  aerothermodynamic  loads  (he.,  mean  and  fluctuating  surface  pressure, 
skin  friction  and  heat  transfer).  Shock  wave  turbulent  boundary  layer  interaction  (SWT- 
BLI)  is  common  in  high  speed  flight,  and  can  significantly  affect  the  aerothermodynamic 
loads.  The  recent  AGARD  Working  Group  18  (WG  18)  study  edited  by  Knight  and  Degrez 
[1]  examined  the  capability  of  general1  Reynolds-averaged  Navier-Stokes  (RANS)  methods 
for  prediction  of  2-D  and  3-D  shock  wave  turbulent  boundary  layer  interactions.  The  study 
comprised  13  configurations  and  37  computations  contributed  by  18  participants.  The  recent 
review  by  Zhcltovodov  [3]  also  examined  the  capability  of  RANS  models  for  prediction  of 
SWTBLI.  Both  studies  concluded  that  general  RANS  methods  accurately  predicted  mean 
surface  pressure  and  heat  transfer,  and  primary  separation  location  for  weak  to  moderate2 
2-D  and  3-D  shock  wave  turbulent  boundary  layer  interactions.  For  strong  interactions,  com¬ 
puted  surface  pressure  and  heat  transfer  distributions,  and  location  of  primary  separation, 
show  significant  disagreement  with  experiment.  In  addition,  the  RANS  calculations  failed  to 
capture  the  high  level  of  unsteadiness  in  the  shock  system  observed  in  the  experiment  and 
the  appearance  of  secondary  separation,  and  no  RANS  calculation  provided  a  prediction  of 
rms  fluctuating  surface  pressure  and  heat  transfer. 

RTO  Working  Group  10  (WG  10)  was  established  in  December  1998  to  examine  three 
specific  technologies  for  propelled  hypersonic  flight:  plug  nozzles  (Subgroup  1),  scramjet 
propulsion  (Subgroup  2)  and  Computational  Fluid  Dynamics  (Subgroup  3).  Six  specific 
topics  were  selected  for  evaluation  by  Subgroup  3,  namely,  boundary  layer  instability  and 
transition,  real  gas  flows,  laminar  hypersonic  viscous-inviscid  interactions,  shock-shock  inter¬ 
actions,  shock  wave  turbulent  boundary  layer  interactions,  and  base  flows  with  and  without 
plume  interactions.  In  each  topic,  specific  experiments  were  identified  for  CFD  validation. 
A  description  of  test  cases  was  published  [4],  The  test  cases  are  also  described  in  [5]. 

This  paper  presents  an  expanded  version  of  the  summary  report  of  the  study  of  shock 
wave  turbulent  boundary  layer  interactions  [6].  Five  separate  configurations  are  considered. 
The  2-D  compression  corner  (Fig.  1)  and  expansion-compression  corner  (Fig.  3)  represent, 
for  example,  deflected  aerodynamic  control  surfaces.  The  2-D  shock  impingement  (Fig.  2) 
typifies  a  shock  boundary  layer  interaction  within  an  inlet,  for  example.  The  3-D  single  fin 
(Fig.  4)  is  typical  of  a  fin-fuselage  juncture.  The  3-D  double  fin  (Fig.  5)  represents  a  sidewall 
compression  inlet.  All  five  configurations  have  been  extensively  studied  and  a  substantial 

1By  general  RANS  method  we  mean  a  RANS  model  which  has  not  been  modified  to  take  into  account 
specific  details  of  the  flowfield  structure.  Significant  progress  has  been  made  in  development  of  engineering 
RANS  methods  which  incorporate  configuration-dependent  flow  physics.  An  example  is  Panaras  [2]. 

2In  this  context,  “weak”  and  “moderate”  refer  to  the  overall  static  pressure  rise. 
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set  of  experimental  data  is  available  for  each  configuration. 

Contributions  of  recent3  numerical  simulations  for  the  five  cases  were  solicited  through 
the  RTO  WG  10  membership  and  the  RTO  WG  10  website  [7].  The  studies  include  Direct 
Numerical  Simulation  (DNS),  Large  Eddy  Simulation  (LES)  and  Reynolds-averaged  Navier 
Stokes  (RANS).  The  computed  results  are  compared  herein  with  the  experimental  data,  and 
the  overall  capabilities  of  the  CFD  methods  are  assessed.  Readers  interested  in  additional 
test  cases  for  SWTBLI  may  consult  Settles  and  Dodson  [8]  and  Haase  et  al  [9]. 


2  2-D  Compression  Corner  and  Shock  Impingement 

The  2-D  supersonic  compression  corner  (Fig.  1)  and  2-D  shock  impingement  (Fig.  2)  are 
standard  CFD  validation  cases  for  shock  wave  turbulent  boundary  layer  interactions.  In 
the  former,  an  equilibrium  supersonic  turbulent  boundary  layer  of  height  6  approaches  a 
compression  corner  defined  by  an  angle  acc.  The  deflection  of  the  flow  by  the  corner  generates 
a  shock  system.  For  sufficiently  large  pressure  rise,  the  boundary  layer  separates  and  a  A- 
shock  forms.  In  the  latter,  a  incident  shock  generated  by  a  wedge  of  angle  aSi  compresses 
the  boundary  layer  leading  to  a  complex  system  of  compression  and  expansion  waves  which 
interact  with  a  reflected  shock. 

The  flowfield  structures  in  the  compression  corner  and  shock  impingement  configurations 
were  investigated  in  the  pioneering  experimental  study  by  Petrov  et  al  [10].  On  the  basis 
of  extensive  experiments  for  laminar  and  turbulent  flows  they  concluded  that  the  2-D  com¬ 
pression  corner  and  2-D  shock  impingement  have  several  similar  features,  and  were  the  first 
to  formulate  the  free  interaction  concept.  In  accordance  with  this  concept,  the  interaction 
in  the  vicinity  of  the  separation  line  is  local  for  developed  separation  zones  and  the  critical 
shock  wave  strength  (the  static  pressure  ratio  across  the  separation  shock  wave)  is  the  same 
for  these  cases  as  well  as  for  the  2-D  shock  wave  turbulent  boundary  layer  interaction  in 
the  vicinity  of  the  shock  wave  which  appears  at  the  exit  in  an  overexpanded  nozzle.  It  does 
not  depend  of  the  strength  of  the  initial  (’’inviscid”)  shock  wave  which  initiates  the  interac¬ 
tion  process  and  depends  only  on  the  freestream  Mach  number  and  the  incoming  boundary 
layer  parameters.  Petrov  et  al  proposed  the  first  empirical  correlation  for  the  critical  shock 
wave  strength  (which  is  close  to  the  specific  plateau  pressure  ratio  in  the  separated  zone) 
and  concluded  that  it  is  approximately  proportional  to  the  freestream  Mach  number  at  high 
Reynolds  numbers.  They  also  proposed  the  first  empirical  correlations  for  the  ratio  of  the 
detachment  distance  of  the  separation  shock  wave  relative  the  surface  trace  of  the  initial 

Simulations  were  selected  from  1997  through  2002,  since  the  review  of  Knight  and  Degrez  [1]  for  AGARD 
WG  18  was  completed  in  early  1997. 
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’’inviscid”  shock  wave  to  the  thickness  of  subsonic  part  of  undisturbed  boundary  layer  for 
the  2-D  compression  corner  and  the  2-D  shock  impingement  configurations.  On  the  basis  of 
these  correlations  they  concluded  that  this  ’’subsonic  region  thickness”  is  the  specific  param¬ 
eter  which  to  a  first  approximation  scales  the  data  (by  a  single  curve  for  every  case)  for  some 
range  of  Mach  and  Reynolds  numbers  as  well  as  the  boundary  layer  thickness  as  a  function 
of  the  flow  inclination  angle  in  the  initial  ’’inviscid”  shock  wave. 

Chapman  et  al  [11]  also  developed  the  free  interaction  concept  through  an  extensive  set  of 
experiments.  Gogish  and  Stepanov  [12]  developed  a  local  theory  of  the  turbulent  boundary 
layer  to  predict  the  critical  shock  wave  strength  and  used  it  to  analyze  the  influence  of  wall 
temperature.  They  successfully  applied  this  theory  to  different  cases  of  turbulent  separation 
including  the  2-D  compression  corner  and  2-D  shock  impingement.  Their  prediction  for  the 
adiabatic  wall  case  is  in  good  agreement  with  the  empirical  correlation  by  Petrov  et  al  [10] 
and  the  free  interaction  theory  of  Chapman  et  al  [11]. 

The  similarity  between  the  flow  properties  for  the  2-D  compression  ramp  and  2-D  shock 
impingement  (as  well  as  for  some  additional  cases),  in  accordance  with  the  free  interaction 
concepts,  were  discussed  and  analyzed  in  the  reviews  by  Green  [13],  Charwat  [14],  Stanewsky 
[15],  Delery  and  Marvin  [16]. 

Similarity  of  the  data  downstream  has  also  been  noted.  Experiments  performed  by  Holden 
[17]  and  Law  [18],  as  well  as  computations  by  Shang  [19]  (Fig.  6)  have  shown  that,  to  a  first 
approximation,  the  surface  pressure  distribution  along  the  entire  length  of  interaction  region 
for  a  compression  corner  of  angle  acc  is  similar  to  a  shock  impingement  generated  by  a  wedge 
of  angle  asi  =  acc/2. 

Despite  such  similarity  in  the  surface  pressure,  differences  in  the  flowfield  structure  and 
curvature  of  the  separated  shear  layer  streamlines  in  these  flows  are  obvious  (Figs.  1  and 
2).  As  described  in  the  monograph  by  Smits  and  Dussuge  [20]  and  in  the  recent  review  by 
Andreopoulos  et  al  [21],  such  differences  are  important  and  can  cause  different  turbulence 
transformation  and  relaxation  processes  in  the  disturbed  boundary  layer  through  the  inter¬ 
action  region.  These  difference  are  manifest  in  the  surface  skin  friction  and  heat  transfer 
distributions  downstream  of  such  interactions. 

2-D  Compression  Corner 

In  his  review,  Delery  [22]  detailed  the  effect  of  the  incoming  boundary  layer  profile  (es¬ 
pecially,  the  incompressible  shape  factor  H{  and  relative  height  ys/S  of  the  sonic  line)  on 
the  size  of  the  separation  region,  and  noted  the  opposite  trends  of  separation  length  (nor¬ 
malized  by  6)  with  Reynolds  number  Res  based  on  the  incoming  boundary  layer  thickness 
6.  In  particular,  for  Res  less  than  approximately  105,  the  separation  length  (normalized  by 
<5)  increases  with  Res  (for  fixed  corner  angle  a  and  wall  temperature  Tu,),  while  the 
opposite  trend  occurs  for  Res  >  105.  Delery  also  summarized  the  effects  of  wall  cooling  and 
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surface  transpiration  on  the  separation  region. 

Zheltovodov  and  Schiilein  [23,  24]  correlated  the  separation  length  Ls  for  a  large  set  of 
experimental  data  of  2-D  compression  corners  according  to 

^  =  f(Res)  (1) 

where 

°  Ml  \Ppl) 

where  P2  is  the  downstream  inviscid  pressure,  and  ppi  is 
according  to  the  empirical  formula  by  Zukoski  [25] 

Ppl  =  Poo  (l-Moo  +  l)  (3) 

Their  correlation,  shown  by  Band  1  in  Fig.  7  from  [24],  encompasses  the  range  3  x  104  < 
Res  <  3  x  106  and  shows  the  change  in  trend  of  separation  length  with  Reynolds  number 
at  Res  ~  105.  The  change  in  behavior  at  Res  ~  105  correlates  with  the  change  [26]  in  the 
exponent  n  of  the  incoming  velocity  prohle 


(2) 


the  plateau  pressure  computed 


u 

uZ 


(4) 


The  variation  of  the  exponent  n  with  Reynolds  number  is  discussed  in  Johnson  and  Bushnell 
[27].  They  showed  that  the  exponent  n  and  velocity  “fullness”  initially  decrease  with  increas¬ 
ing  Reynolds  number,  but  that  the  opposite  behavior  occurs  at  higher  Reynolds  number. 
Increased  values  of  n  implies  a  higher  resistance  to  separation  (see  also  Delery  and  Marvin 
[16]). 

Settles  [28,  29]  developed  an  expression  for  the  variation  of  the  upstream  influence  length 
Lj,  defined  as  the  distance  between  the  initial  pressure  rise  and  the  corner,  with  Reynolds 
number  Res  and  a  at  Moo  ~  3, 

=  0.9  e0'23"  (5) 

o 

The  correlation  was  verified  for  105  <  Res  <  10'  and  10°  <  a  <  26°. 

Adamson  and  Messiter  [30]  described  the  application  of  Triple  Deck  Theory,  originally  de¬ 
veloped  for  separation  of  a  laminar  compressible  boundary  layer  by  Stewartson  and  Williams 
[31]  and  Neiland  [32],  to  turbulent  boundary  layer  separation  associated  with  shock  impinge¬ 
ment.  Agrawal  and  Messiter  [33]  analyzed  an  unseparated  2-D  compression  corner  using 
asymptotic  methods. 
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Zheltovodov  [34]  described  the  principal  elements  of  the  flowfield  structure  and  turbulence 
transformation  in  the  compression  corner  (Fig.  8)  on  the  basis  of  a  series  of  experimental 
studies  [35,  36,  37,  38,  39,  40,  41]  which  included  optical  (spark  shadowgraph  and  Schlieren) 
flowfield  visualization  and  surface  flow  pattern  visualization,  surface  pressure,  skin  friction 
and  heat  transfer  measurements  as  well  as  pitot  and  static  pressure  survey,  velocity  profiles, 
hot  wire  measurements  of  the  mass  flow,  velocity,  density  and  total  temperature  fluctuations. 
These  elements  include  the  amplification  of  the  turbulence  by  shock  waves  in  the  boundary 
layer  (1)  and  external  flow  (2),  suppression  of  turbulence  by  the  expansion  fan  (3),  formation 
of  a  new  layer  in  the  near  wall  part  of  the  attaching  flow  (4),  formation  of  Taylor-Gortler 
vortices  (5),  and  reverse  transition  (6)  in  the  separation  region  due  to  the  favorable  pressure 
gradient  and  decreasing  of  the  local  Reynolds  number  in  the  reverse  flow  caused  by  the 
decrease  in  the  velocity  in  the  direction  of  the  separation  point. 

The  unsteadiness  of  the  shock  system  and  separation  flow  in  the  compression  corner  is 
an  additional  important  aspect.  Detailed  reviews  have  been  prepared  by  Dolling  [42,  43], 
Smits  and  Dussauge  [20]  and  Adams  [44],  This  phenomena  have  been  analyzed  by  several 
investigators.  For  example,  the  shock  wave  deformations  cased  by  the  unsteadiness  have  been 
observed  using  the  optical  (spark  Schlieren)  visualization  by  Zheltovodov  et  al  [36,  45]  and 
using  hot  wires  by  Zheltovodov  and  Yakovlev  [35,  46].  It  was  shown  that  the  appearance  of 
separation  and  accompanying  large  scale  eddies  in  the  separated  shear  layer  for  a  =  25°  and 
45°  at  Mach  2.9  stimulates  a  very  high  level  of  fluctuations  of  mass  flow,  velocity  and  total 
temperature  (compared  to  the  unseparated  interaction  at  a  =  8°)  in  both  the  separated  and 
attached  boundary  layers  and  in  the  external  flow.  For  example,  the  mass  flow  fluctuations 
in  the  external  flow  above  the  boundary  layer  edge  at  a  =  25°  and  45°  were  amplified  by 
a  factor  of  11  to  12,  while  the  corresponding  amplification  factor  was  approximately  2  for 
a  =  8°. 

Zheltovodov  et  al  [34,  39,  40,  47]  and  Borisov  et  al  [48]  concluded  on  the  basis  of  experi¬ 
mental  studies  and  theoretical  analysis  that  the  amplification  of  turbulence  is  the  important 
factor  which  stimulates  the  increase  in  the  surface  heat  transfer  level  in  the  interaction  region 
as  well  as  downstream  in  the  attached  flow.  The  DNS  of  Adams  [44]  at  Mach  3  and  a  =  18° 
revealed  compression  waves  emanating  from  the  separated  boundary  layer  edge  which  are  in 
qualitative  agreement  with  the  Schlieren  visualization  of  Zheltovodov  et  al  at  Mach  2.9  and 
a  =  25°.  Both  computations  and  experiments  therefore  indicate  that  moving  compression 
waves  are  formed  at  the  upstream  edges  of  the  large-scale  eddies  which  appear  in  the  process 
of  oscillations  of  the  separation  shock. 

The  unsteadiness  of  the  shock  system  in  the  compression  corner  has  been  studied  in  detail 
by  Dolling  and  Murphy  [49]  and  Dolling  and  Or  [50]  who  observed  separation  shock  motion 
on  the  scale  of  5  at  frequencies  more  than  an  order  of  magnitude  below  the  frequency  Uoo/5 
of  the  energy  containing  eddies  of  the  incoming  boundary  layer.  Andreopoulos  and  Muck 
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[51]  observed  that  the  frequency  of  the  separation  shock  unsteadiness  is  approximately  the 
same  as  the  bursting  frequency  of  the  upstream  boundary  layer.  Erengil  and  Dolling  [52,  53] 
identified  both  small-  and  large-scale  fluctuations  of  the  separation  shock  motion.  Dolling 
[54]  reviewed  the  literature  on  shock  unsteadiness  through  the  early  1990s;  see  also  Dolling 
[42,  43], 

Bibko  et  al  [55,  56]  performed  detailed  experimental  studies  of  surface  pressure  fluctuation 
distributions  as  well  as  their  spectra  upstream  of  compression  corners  for  different  inclination 
angles  a  =  10°  to  90°  in  the  Mach  number  range  from  2  to  4.  They  concluded  that  the 
influence  of  the  incoming  turbulent  boundary  layer  fluctuations  on  the  fluctuations  of  the 
shock  wave  induced  by  separation  (the  ’’separation  shock”)  is  very  small.  According  to  their 
investigations,  the  main  reason  for  the  shock  wave  fluctuations  is  the  disturbances  which 
appear  in  the  developed  (large  scale)  separation  zone  as  well  as  in  the  attachment  regions. 
On  a  basis  of  analysis  of  the  experiments  by  Bibko  et  al,  Zheltovodov  [3]  concluded  that  the 
stabilization  of  the  surface  pressure  fluctuation  spectra  in  the  separation  zone  as  well  as  in 
the  vicinity  of  the  separation  shock  with  increasing  a  is  attributable  to  the  inviscid  shock 
detachment  when  the  maximum  separation  zone  is  formed. 

Glotov  [57]  proposed  a  physical  interpretation  of  the  shock  wave  unsteadiness  variation 
in  the  vicinity  of  axisymmetric  compression  corners  at  different  stages  of  the  separation 
development  in  the  process  of  its  interaction  with  the  large  eddies.  Using  simultaneous 
high  speed  shadowgraphy  and  surface  pressure  fluctuation  measurements  in  the  vicinity 
of  the  corner  apex,  he  concluded  that  for  small  scale  separation  the  large  eddies  in  the 
external  part  of  the  boundary  layer  periodically  supplant  the  separation  (reverse  flow)  zone 
and  prevent  the  penetration  of  disturbances  from  downstream  to  upstream.  He  concluded 
that  this  process  stimulates  large  deformations  of  the  shock  wave  similar  to  the  interaction 
of  a  shock  wave  with  a  vortex  ring.  For  large  scale  separation,  such  eddies  are  displaced 
downstream  of  the  separation  shock  wave  above  the  surface  together  with  the  separated  shear 
layer,  and  the  disturbances  penetrating  upstream  from  the  attachment  region  play  the  main 
role  in  stimulating  low  frequency  shock  wave  fluctuations,  with  high  frequency  fluctuations 
stimulated  by  the  former  mechanism  superimposed  upon  these.  This  interpretation  is  in 
agreement  with  conclusions  by  Bibko  et  al. 

Wu  and  Miles  [58]  visualized  the  interaction  of  the  large  eddies  in  the  incoming  bound¬ 
ary  layer  with  the  separation  shock  using  a  high  speed  (MHz)  pulse-burst  laser  with  a 
charge-coupled-device  framing  camera.  They  observed  shock  motion  on  the  scale  of  6  at 
the  frequency  of  the  incoming  eddies  Uoo/6  and  compression  of  the  eddies  in  the  streamwise 
direction  by  the  shock. 

Zheltovodov  [3,  47]  and  Dolling  [42]  summarized  the  status  of  RANS  simulations  of  the 
compression  corner  flow.  They  observed  that  no  RANS  model  was  capable  of  accurate 
prediction  of  all  aspects  of  the  flowfield.  In  particular,  the  RANS  calculations  did  not 
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display  the  unsteadiness  of  the  separation  shock  which  was  concluded  to  be  a  major  defect 
in  the  simulation  and  the  cause  of  the  lack  of  agreement  with  the  experimental  mean  surface 
pressure. 

Yan  et  al  [59]  compiled  a  summary  of  RANS  predictions  of  separation  length  for  various 
2-D  compression  corner  configurations  [60,  61,  62,  63,  64,  65].  In  general  the  simulations 
(Fig.  9,  see  also  Table  1)  do  not  accurately  predict  the  separation  length  in  comparison  with 
the  experimental  correlation  (1). 

Ad  hoc  modifications  to  RANS  models  have  achieved  improved  prediction  of  the  2-D 
compression  corner;  however,  such  modifications  are  not  general  in  nature  and  therefore 
cannot  be  readily  extended  to  other  SWTBLI  configurations.  Examples  are  Borisov  et  al 
[65]  and  Bedarev  et  al  [66]  who  demonstrated  the  capability  for  the  k  —  u  model  to  predict 
relaminarization  in  the  reversed  flow  region  through  an  ad  hoc  modification  of  uj.  These 
computations  indicate  the  necessity  to  model  this  phenomenon  in  order  to  achieve  better 
prediction  of  the  surface  skin  friction  and  heat  transfer  in  the  separation  zone. 

2-D  Shock  Impingement 

The  studies  of  2-D  shock  impingement  are  not  as  numerous  as  the  2-D  compression 
corner.  Several  of  these  studies  have  been  reviewed  by  Green  [13],  Charwat  [14],  Stanewsky 
[15],  Delery  and  Marvin  [16]. 

Kuehn  [67]  performed  a  series  of  experiments  for  2-D  shock  impingement  at  Mach  2.6  to 
2.85  and  proposed  criteria  for  incipient  separation  based  on  the  appearance  of  an  inflection 
point  in  the  mean  surface  pressure  distribution.  He  also  observed  unsteadiness  in  the  shock 
structure  for  separated  interactions. 

Reda  and  Murphy  [68]  measured  surface  pressure  and  pitot/static  pressure  profiles  for  a 
2-D  shock  impingement  at  Mach  2.9  for  a  =  5°  to  13°.  Significant  3-D  effects  were  observed 
in  both  surface  pressure  and  oil  him  visualization  due  to  test  section  side  wall  and  corner 
boundary  layer  effects.  Reda  and  Murphy  [69]  examined  the  effect  of  the  sidewall  boundary 
layers  on  incipient  separation  and  interaction  size. 

Law  [70]  measured  surface  pressure  and  howheld  density  (using  a  holographic  interfer¬ 
ometer  technique)  at  Mach  2.96  for  a  =  8°  to  11°.  The  holographic  images,  however, 
incorporated  significant  3-D  effects  within  and  downstream  of  the  interaction  region  due  to 
corner  effects. 

Kussoy  and  Horstman  [71]  performed  experiments  for  an  axisymmetric  shock  impinge¬ 
ment  at  Mach  7  for  a  =  7.5°  and  15°.  Measurements  include  surface  pressure,  temperature 
and  skin  friction,  and  profiles  of  pitot  pressure,  static  pressure  and  total  temperature  at 
several  locations.  Kussoy  and  Horstman  [72]  measured  surface  pressure  and  heat  transfer 
for  a  2-D  shock  impingement  interaction  at  Mach  8.2  for  a  =  5°  to  10°. 

Zakharov  and  Kutuzova  [73]  developed  approximate  (semi-empirical)  methods  to  predict 
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the  upstream  influence  line  position  (the  line  where  the  surface  pressure  starts  to  increase) 
regarding  the  trace  of  "inviscid”  shock  for  its  interaction  with  laminar,  transitional  and 
turbulent  boundary  layers. 

Brazhko  [74]  performed  optical  visualization  of  Goertler-type  vortices  using  the  laser 
sheet  technique  and  showed  that  they  resulted  in  periodical  variation  of  the  heat  transfer  at 
Mach  5  and  6  for  both  laminar  and  turbulent  regimes  of  separated  flows. 

The  studies  by  Rose  and  Johnson  [75],  Mikulla  and  Horstman  [76],  Modarress  and  John¬ 
son  [77],  and  Meyer  et  al  [78]  provide  detailed  information  regarding  the  mean  flowfield 
and  turbulence  measurements.  They  demonstrate  the  increase  in  fluctuations  of  mass  flow 
and  velocity,  turbulence  intensity,  Reynolds  shear  stress  and  turbulent  kinetic  energy  down¬ 
stream  of  the  interaction.  The  strength  of  the  impingement  shock  has  a  significant  effect 
on  the  downstream  properties.  Most  of  the  turbulence  properties  also  increased  with  shock 
strength  near  the  interaction  in  unseparated  cases  with  the  exception  of  the  cross-stream 
velocity  fluctuations.  It  was  also  noted  that  the  mean  flow  and  turbulent  properties  as  well 
as  relaxation  of  the  disturbed  flow  aft  of  the  separation  region  were  noticeably  different 
than  observed  for  unseparated  cases.  A  strong  influence  of  the  separation  flow  unsteadiness 
on  the  turbulence  energy  was  revealed.  The  increase  of  turbulence  intensity  and  Reynolds 
shear  stress  downstream  of  the  interaction  in  such  situations  was  small  compared  to  the 
unseparated  case. 

The  complex  influence  of  the  flowfield  unsteadiness  may  explain  the  differences  between 
RANS  computations  and  experimental  data  when  separation  occurs.  For  example,  recent 
computations  by  Fedorova  et  al  [79]  using  the  k  —  u  turbulence  model  clearly  demonstrate 
increasing  differences  between  computed  and  experimental  velocity  profiles,  surface  skin 
friction  and  heat  transfer  with  stronger  impinging  shock  waves.  Linbland  et  al  [80]  and 
Rizzetta  [81]  showed  significant  differences  in  the  prediction  of  flow  properties  using  different 
turbulence  models.  These  conclusions  are  supported  by  the  measurements  of  fluctuating 
heat  transfer  performed  by  Hayashi  et  al  [82]  for  both  unseparated  and  separated  flows,  as 
well  as  the  subsequent  analysis  of  this  data  by  Brusniak  and  Dolling  [83].  Their  analysis 
shows  that  the  fluctuating  heat  transfer  signals  from  a  region  of  shock  motion  bear  a  close 
resemblance  to  the  character  of  fluctuating  pressure  signals  taken  from  the  same  region. 
Analysis  of  the  data  for  different  cases  including  the  2-D  compression  corner  and  the  2-D 
shock  impingement  demonstrates  the  important  role  of  shock  unsteadiness  and  indicates 
directions  for  development  of  adequate  engineering  and  CFD  methods  for  prediction. 

Holden  and  Chadwick  [84]  conducted  an  experimental  study  of  the  effects  of  tangential 
film  and  transpiration  cooling  on  surface  heat  transfer  in  a  2-D  shock  impingement  at  Mach  6 
and  8.  They  concluded  that  film  cooling  was  ineffective  for  all  but  weak  shock  interactions 
due  to  the  dispersal  of  the  film  by  the  shock  boundary  layer  interaction.  Transpiration 
cooling  was  found  to  be  effective. 
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In  the  following  sections,  DNS  and  LES  results  are  presented  for  several  configurations 
of  the  2-D  compression  corner  and  one  configuration  of  the  2-D  shock  impingement.  Mean 
quantities  represent  combined  time  and  spanwise  averages. 


2.1  Case  1:  Mach  3,  a  =  8° 


Urbin  et  al  [85]  and  Kannepalli  et  al  [86]  simulated  a  Mach  3,  a  =  8°  compression  corner 
at  Res  =  2  x  104  and  adiabatic  wall  conditions.  Experimental  data  of  Zheltovodov  et  al 
[35,  36,  37]  is  available  for  a  =  8°  at  Reynolds  numbers  Res  =  7.5  x  104  and  1.1  x  105. 
The  flow  conditions  are  summarized  in  Table  2.  The  inviscid  pressure  ratio  across  the  shock 
P2/P1  =  1-8.  In  accordance  with  the  boundaries  of  the  different  stages  of  separation  develop¬ 
ment  in  compression  ramps  as  defined  by  Zheltovodov  [3],  this  configuration  corresponds  to 
an  incipient  intermittent  separation.  Although  the  mean  flowheld  is  not  separated,  instan¬ 
taneous  reversed  flow  can  exist  for  a  short  time  in  the  vicinity  of  the  shock  wave  and  near 
the  wall  due  to  the  unsteadiness  of  the  interaction  process. 

Details  of  the  computations  are  presented  in  Table  3  where  SA  and  TA  indicate  the 
spatial  and  temporal  accuracy.  Urbin  et  al  used  an  unstructured  grid  LES  algorithm  with  a 
static  Smagorinsky  model.  The  inflow  conditions  were  obtained  from  a  separate  computation 
of  a  flat  plate  adiabatic  boundary  layer  at  Mach  3.  The  inflow  velocity  profile  satisfied  the 
compressible  Law  of  the  Wall  [87].  The  inflow  skin  friction  coefficient  agreed  with  the 
empirical  value  obtained  from  the  compressible  Law  of  the  Wall  and  Wake  to  within  6%, 
and  the  inflow  adiabatic  wall  temperature  was  within  3%  of  the  theoretical  value  obtained 
from  the  empirical  formula 


T  =  T 

a. in  rx 


1  + 


(7-1) 


PnMi 


(6) 


where  Prtm  =  0.89  is  the  mean  turbulent  Prandtl  number.  Kannepalli  et  al  used  the  MILES 
[88]  (Monotone  Integrated  Large  Eddy  Simulation)  approach  wherein  the  inherent  dissipation 
in  the  numerical  algorithm  constitutes  the  SGS  model. 

An  instantaneous  image  of  the  shock  wave  (as  visualized  by  the  isosurface  p/p0 0  =  1.25) 
and  the  streamwise  velocity  at  y+  =  10  is  shown  in  Fig.  10  from  Urbin  et  al.  Regions  of 
negative  velocity  (instantaneous  separation)  are  visible  in  agreement  with  the  conclusion  of 
Zheltovodov  [3]  regarding  an  incipient  intermittent  separation  regime  at  such  conditions. 

The  mean  wall  pressure  distribution  is  shown  in  Fig.  11  for  Urbin  et  al.  The  experimental 
profiles  for  the  two  different  Reynolds  numbers  are  essentially  identical.  The  computed  profile 
shows  close  agreement  with  experiment. 

Additional  results  by  Urbin  et  al  are  presented  in  Figs.  12  to  15.  The  mean  velocity 
profiles  upstream  and  downstream  of  the  corner  are  displayed  in  Fig.  12.  The  experimental 
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data  is  for  Res  =  1.1  x  105.  The  upstream  profile  is  3.35  ahead  of  the  corner  and  in  the 
undisturbed  flat  plate  boundary  layer.  The  downstream  profile  is  3.55  from  the  corner  and 
located  close  to  the  position  where  the  wall  pressure  reaches  its  asymptotic  value.  The 
computed  results  are  in  close  agreement  with  experiment.  The  mean  static  temperature 
profiles  at  the  same  locations  are  shown  in  Fig.  13.  Upstream  of  the  interaction  the  computed 
and  experimental  profiles  are  in  close  agreement.  Downstream  of  the  interaction,  the  profiles 
agree  except  within  the  region  0.05  <  y/5  <  0.25  where  the  experimental  value  exceeds  the 
prediction. 

The  Reynolds  streamwise  and  shear  stresses  are  displayed  in  Figs.  14  and  15,  respectively, 
where  the  stresses  have  been  normalized  by  the  upstream  wall  shear  stress  rw.  The  experi¬ 
mental  data  of  Muck  et  al  [89]  is  shown.  The  experimental  Reynolds  number  Res  =  1.6  x  106 
is  a  factor  of  eighty  greater  than  the  LES.  However,  scaling  of  the  Reynolds  shear  stress  in  a 
flat  plate  boundary  layer  by  the  wall  shear  stress  rw  has  been  shown  to  collapse  the  exper¬ 
imental  data  for  a  range  of  Mach  and  Reynolds  numbers  [20].  Thus,  the  comparison  of  the 
LES  and  experiment  is  justified  if  only  for  qualitative  comparison.  The  effect  of  the  shock 
boundary  layer  interaction  is  to  increase  the  level  of  Reynolds  normal  pu"u"  and  shear  pu"v" 
stresses.  The  predicted  values  are  in  good  agreement  with  experiment  except  for  y/5  <  0.1. 

In  summary,  the  following  conclusions  can  be  made: 

•  LES  of  the  Mach  3,  a  =  8°  compression  corner  shows  good  agreement  with  the  experi¬ 
mental  mean  surface  pressure  distribution  and  qualitatively  correct  modification  of  the 
mean  velocity  and  turbulence  profiles  by  the  shock  wave 

•  The  mean  flowfield  is  unseparated  (and  consequently,  this  case  does  not  represent  a 
particularly  difficult  case  for  LES),  but  the  instantaneous  flowfield  indicates  an  incipient 
intermittent  separation  regime  in  accordance  with  the  generalization  of  experimental 
data  and  analysis  by  Zheltovodov  [3] 

2.2  Case  2:  Mach  3,  a  =  18° 

Adams  [44,  90,  91],  Rizzetta  et  al  [92]  and  Stolz  et  al  [93]  simulated  a  Mach  3,  a  =  18° 
compression  corner  at  Res  =  2.1  x  104  and  adiabatic  wall  conditions.  The  flow  conditions 
are  summarized  in  Table  4.  No  experimental  data  is  available  for  this  configuration.  However, 
in  accordance  with  the  generalizations  by  Zheltovodov  [3]  for  the  boundaries  of  the  different 
stages  of  separation  development,  mean  separation  occurs  for  this  configuration.  The  inviscid 
pressure  ratio  across  the  shock  is  p2/Pi  =  3.37. 

Details  of  the  computations  are  presented  in  Table  5.  Adams  performed  a  DNS.  Rizzetta 
et  al  used  three  different  approaches,  namely,  a  DNS,  an  LES  using  a  static  Smagorin- 
sky  model  for  the  SGS  stresses  (LES-S),  and  an  LES  using  a  dynamic  (Germano  et  al 
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[94])  Smagorinksy  model  (LES-D).  Stolz  et  al  used  an  Approximate  Deconvolution  Method 
(ADM).  The  inflow  conditions  were  obtained  in  each  instance  from  a  separate  computation 
of  a  flat  plate  adiabatic  boundary  layer  at  Mach  3  (Fig.  16). 

A  comparison  of  the  simulations  of  Adams  and  Rizzetta  et  al  are  presented  in  Figs.  17  to 
24.  Although  the  Mach  number  and  Reynolds  number  were  nominally  identical,  the  inflow 
profiles  differ  significantly.  Figs.  17  and  18  display  the  mean  and  rms  fluctuating  streamwise 
velocity  at  the  inflow  boundary.  Major  differences  are  evident,  and  have  a  significant  effect 
on  the  computed  shock  boundary  layer  interaction  as  discussed  below. 

The  computed  mean  surface  pressure  coefficient  cv  =  (p— Poc) / 1 PocU<L  is  shown  in  Fig.  19. 
The  size  of  the  interaction  region,  defined  by  the  location  of  the  initial  rise  in  cp  and  inter¬ 
section  of  cp  with  the  downstream  inviscid  value  (cp  =  0.376),  is  substantially  larger  in  the 
simulation  of  Adams  compared  to  Rizzetta.  This  is  attributable  to  the  differences  in  the 
inflow  velocity  profile  (Fig.  17),  as  the  momentum  flux  pu  is  lower  near  the  wall  in  Adams’ 
profile  (see  also  Fig.  22).  The  DNS  and  LES  profiles  of  Rizzetta  are  virtually  identical, 
implying  the  SGS  model  has  negligible  effect. 

The  computed  mean  skin  friction  coefficient  Cf  =  Twl\p00U‘l0  is  shown  in  Fig.  20.  The 
upstream  skin  friction  level  is  virtually  identical  for  both  Adams  and  Rizzetta  (see  Fig.  17 
which  indicates  the  velocity  profiles  are  essentially  the  same  within  the  viscous  sublayer). 
The  separation  length  predicted  by  Adams  is  approximately  130%  larger  than  Rizzetta.  This 
is  attributable  to  the  lower  momentum  flux  in  the  incoming  boundary  layer  (Fig.  17).  The 
DNS  and  LES  profiles  of  Rizzetta  are  virtually  identical,  again  indicating  the  SGS  model 
has  negligible  effect. 

The  mean  velocity  and  static  temperature  are  displayed  in  Figs.  21  and  22,  respectively. 
Station  Nos.  1,  2  and  3  are  located  upstream  of  the  interaction  region,  at  the  corner,  and 
downstream  of  the  corner,  respectively.  The  mean  velocity  profiles  are  displayed  using  the 
Van  Driest  [95]  transformation.  All  profiles  exhibit  agreement  with  the  Law  of  the  Wall 
logarithmic  profile  in  the  region  20  <  y+  <  100  (Station  No.  1).  Downstream  of  the  corner, 
the  simulations  by  Rizzetta  show  recovery  of  the  near  wall  region  to  the  equilibrium  Law 
of  the  Wall,  while  the  simulation  by  Adams  does  not.  Significant  differences  in  the  mean 
temperature  between  Adams  and  Rizzetta  are  evident  both  upstream,  at  and  downstream 
of  the  corner.  These  are  associated  with  differences  in  the  mean  velocity,  since  the  total 
temperature  is  approximately  constant  across  the  boundary  layer.  For  both  the  velocity  and 
temperature  profiles,  the  DNS  and  LES  results  of  Rizzetta  are  virtually  identical. 

Spanwise  turbulent  kinetic  energy  spectra  at  two  separate  positions  within  the  boundary 
layer  and  three  streamwise  locations  are  shown  in  Figs.  23  and  24  where  the  location  is 
defined  in  terms  of  the  inflow  boundary  layer  condition.  All  curves  have  been  normalized 
to  unity  at  k$  =  1.  The  results  of  Rizzetta  are  consistently  higher  than  Adams,  implying  a 
higher  level  of  turbulent  kinetic  energy.  This  result  is  confirmed  by  comparison  of  profiles  for 
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the  streamwise,  normal  vertical  and  spanwise  rms  fluctuating  velocities  [92] .  The  profiles  of 
Rizzetta  display  the  emergence  of  an  inertial  subrange  (E(k)  oc  £v5/3).  The  inertial  subrange 
properly  occurs  at  wavenumbers  above  the  peak  in  the  energy  spectrum  [96].  Adams’  results 
do  not  display  an  inertial  subrange.  The  spectra  of  Rizzetta  at  high  wavenumber  show  higher 
energy  for  the  LES  than  DNS  calculations.  This  is  an  unexplained  result,  since  it  would  be 
expected  that  the  SGS  model  would  dissipate  more  energy  at  higher  wavenumber  than  the 
DNS. 

A  comparison  of  the  results  of  Adams  and  Stolz  et  al  show  close  agreement.  An  example 
is  presented  in  Fig.  25  which  displays  the  rms  of  the  Reynolds  fluctuations  of  the  streamwise 
momentum,  density,  streamwise  velocity  and  static  temperature.  The  maximum  difference 
between  the  DNS  and  LES  is  less  than  5%  everywhere.  The  LES,  however,  required  sub¬ 
stantially  less  computer  resources. 

The  computed  separation  length,  normalized  using  (2),  is  shown  in  Fig.  26.  The  computed 
separation  lengths  of  Rizzetta  are  below  the  linear  extrapolation  of  the  lower  bound  of  the 
experimental  band,  while  the  computed  separation  length  of  Adams  and  Stolz  et  al  lies  above 
the  upper  bound. 

In  summary,  the  following  conclusions  may  be  made: 

•  DNS  of  shock  wave  turbulent  boundary  layer  interactions  is  sensitive  to  the  inflow 
profiles 

•  The  DNS  and  LES  of  Rizzetta  are  virtually  identical  for  all  profiles.  This  implies  that 
the  SGS  model  has  a  negligible  effect. 

•  The  computed  mean  separation  lengths  are  outside  the  linear  extrapolation  of  the 
experimental  data 

•  The  DNS  of  Adams  and  LES  (ADM)  of  Stolz  et  al  are  in  close  agreement 

2.3  Case  3:  Mach  3,  a  =  24° 

Rizzetta  and  Visbal  [97]  and  Kannepalli  et  al  [86]  simulated  a  Mach  3,  a  =  24°  compression 
corner  at  Res  =  2  x  104  and  adiabatic  wall  conditions4.  Experimental  data  is  available  from 
Settles  et  al  [98],  Dolling  and  Murphy  [49]  and  Srnits  and  Muck  [99]  at  a  =  24°  and  from 
Zhcltovodov  et  al  [3,  35,  36,  37,  100]  at  a  =  25°  as  indicated  in  Table  6.  The  experimental 
data  of  Zhcltovodov  et  al  at  a  =  25°  are  referenced  here  since  Kannepalli  et  al  compared 
their  results  at  a  =  24°  with  this  experimental  dataset.  The  Reynolds  numbers  for  the 
experiments  at  a  =  24°  are  significantly  higher  than  the  simulation  of  Rizzetta  and  Visbal 

4Rizzetta  and  Visbal  [97]  also  simulated  a  =  8°,  16°  and  20°  at  Res  =  2x  104. 
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by  a  factor  of  40  to  77,  and  therefore  the  comparison  between  this  simulation  and  experiment 
is  qualitative.  The  Reynolds  number  for  the  experiments  at  a  =  25°  are  a  factor  of  3.1  larger 
than  the  simulation  of  Kannepalli  et  al.  The  inviscid  pressure  ratio  across  the  shock  for  the 
a  =  24°  corner  is  P2/P1  =  4.68. 

Details  of  the  computations5  are  presented  in  Table  7.  A  structured  grid  LES  algorithm 
using  a  dynamic  Smagorinsky  model  was  utilized  by  Rizzetta  and  Visbal,  while  the  MILES 
methodology  was  utilized  by  Kannepalli  et  al. 

The  results  of  Rizzetta  and  Visbal  are  presented  in  Figs.  26  to  35.  The  effect  of  the 
difference  in  Reynolds  number  is  evident  in  the  inflow  profiles  of  streamwise  velocity  (Fig.  27), 
streamwise  Reynolds  stress  (Fig.  28)  and  Reynolds  shear  stress  (Fig.  29).  The  computed 
inflow  mean  velocity  is  lower  than  the  experiment  near  the  wall  as  expected. 

The  mean  wall  pressure  is  shown  in  Fig.  30.  The  upstream  influence,  plateau  and  re¬ 
covery  regions  are  accurately  predicted,  notwithstanding  the  difference  in  Reynolds  number 
between  the  simulation  and  experiment.  The  standard  deviation  ap  of  the  surface  pressure 
fluctuations,  normalized  by  the  local  mean  wall  pressure  pw,  is  displayed  in  Fig.  31.  The  ex¬ 
perimental  [49]  peak  value  is  associated  with  the  fluctuating  motion  of  the  separation  shock. 
The  simulation  significantly  underestimates  the  peak  value,  and  overestimates  the  level  of 
fluctuations  downstream  of  the  corner.  Similarly,  Fig.  32  indicates  that  the  simulation  un¬ 
derestimates  the  peak  value  of  skewness  S  =  p/3/cr3  near  separation.  These  discrepancies 
can  be  understood  by  considering  the  time  history  of  the  surface  pressure  fluctuations  in  the 
experiment  and  simulation.  Fig.  33  shows  a  time  history  of  the  surface  pressure  at  midspan 
and  x/5  =  —2.1  which  is  near  the  location  of  the  peak  values  of  the  standard  deviation 
(Fig.  31)  and  skewness  (Fig.  32).  The  experiment  displays  a  high  intermittency  of  the  pres¬ 
sure  signal  associated  with  low  frequency  motion  of  the  separation  shock.  The  simulation 
displays  no  low  frequency  component. 

The  mean  skin  friction  coefficient  is  shown  in  Fig.  34.  The  position  of  separation  and 
reattachment  (as  measured  from  the  corner  at  x  —  0)  are  significantly  overestimated,  and 
the  recovery  of  the  skin  friction  downstream  of  attachment  substantially  underpredicts  the 
experiment.  The  computed  separation  length,  normalized  using  (2),  is  shown  in  Fig.  26.  The 
computed  separation  length  is  at  the  lower  limit  of  a  linear  extrapolation  of  the  experimental 
band.  Velocity  profiles  at  two  locations  downstream  of  the  corner  at  x/5  —  1.33  and  2.65  are 
displayed  in  Fig.  35.  The  underestimate  in  the  velocity  in  the  near  wall  region  is  consistent 
with  the  comparison  of  computed  and  experimental  upstream  profiles  (Fig.  27). 

The  results  of  Kannepalli  et  al  are  presented  in  Figs.  36  and  37.  The  mean  surface 
pressure  (Fig.  36)  shows  good  agreement  with  the  experimental  data  of  Zhcltovodov  et  al 

5A  coarse  grid  simulation  (211  x  76  x  81)  was  also  performed  by  Rizzetta  and  Visbal  which  demonstrated 
that  the  fine  grid  solution  for  the  surface  pressure  was  grid  converged. 
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including  the  presence  of  the  plateau  region.  The  computed  results  of  Yan  et  al  [59]  at 
a  =  25°,  discussed  in  the  next  section,  are  shown  for  comparison.  The  mean  skin  friction 
(Fig.  37)  also  shows  good  agreement  with  experiment,  and  in  particular  is  in  good  agreement 
with  the  mean  separation  and  reattachment  locations.  The  computed  separation  length, 
normalized  using  (2),  is  shown  in  Fig.  26  and  is  within  the  limits  of  a  linear  extrapolation 
of  the  experimental  band. 

In  summary,  the  following  conclusions  can  be  made: 

•  The  mean  surface  pressure  predicted  by  Rizzetta  and  Visbal  shows  good  agreement 
with  experiment.  This  may  be  fortuitous,  however,  in  view  of  the  significant  differ¬ 
ence  in  Reynolds  number  between  the  simulation  and  experiment,  and  the  substantial 
differences  in  surface  pressure  fluctuation  statistics  in  the  vicinity  of  separation.  The 
computed  separation  length  is  at  the  lower  limit  of  the  experimental  data  correlation. 
The  near  wall  velocity  profiles  downstream  of  the  corner  are  significantly  underpre¬ 
dicted,  likely  due  to  the  difference  in  inflow  profiles  associated  with  the  disparity  in 
Reynolds  numbers. 

•  The  mean  surface  pressure  and  skin  friction  predicted  by  Kannepalli  et  al  show  good 
agreement  with  the  experiment.  The  Reynolds  numbers  of  the  simulation  and  experi¬ 
ment  are  much  closer  than  for  Rizzetta  and  Visbal. 

2.4  Case  4:  Mach  3,  a  =  25° 

Yan  et  al  [59]  simulated  a  Mach  2.88,  a  =  25°  compression  corner  at  Res  =  2  x  104  and  adia¬ 
batic  wall  conditions.  Experimental  data  of  Zhcltovodov  et  al  [3,  35,  36,  37,  100]  is  available 
for  a  =  25°  at  Reynolds  number  Res  =  6.4  x  104.  The  flow  conditions  are  summarized  in 
Table  8.  The  inviscid  pressure  ratio  across  the  shock  p2/fh  =  4.7.  The  mean  flowfield  is 
separated  at  the  corner. 

Details  of  the  computation  are  presented  in  Table  9.  An  unstructured  grid  LES  algorithm 
using  the  MILES  model  [88]  was  utilized.  The  inflow  velocity  profile  satisfied  the  Law  of  the 
Wall.  The  inflow  skin  friction  coefficient  agreed  with  the  empirical  value  obtained  from  the 
compressible  Law  of  the  Wall  and  Wake  to  within  6%.  The  inflow  adiabatic  wall  temperature 
was  within  3%  of  the  theoretical  value  obtained  from  (6). 

The  mean  wall  pressure  is  shown  in  Fig.  38.  The  upstream  influence  (i.e.,  the  location  of 
the  initial  pressure  rise)  is  accurately  predicted;  however,  the  computed  pressure  distribution 
does  not  show  the  distinct  plateau  evident  in  the  experiment,  and  recovers  to  the  downstream 
inviscid  pressure  more  rapidly  than  the  experiment. 

The  mean  skin  friction  coefficient  is  shown  in  Fig.  39.  The  initial  drop  in  skin  friction 
associated  with  the  adverse  pressure  gradient  is  reasonably  predicted.  The  predicted  recovery 
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of  the  skin  friction  downstream  of  reattachment  is  more  rapid  than  in  the  experiment,  similar 
to  the  results  for  the  surface  pressure.  The  computed  mean  separation  length  is  shown  in 
Fig.  26.  The  value  lies  within  a  linear  extrapolation  of  the  band  of  experimental  data. 

Fig.  40  displays  the  instantaneous  shock  wave  structure  defined  by  isosurfaces  of  static 
pressure  at  p/poo  =  1-4  and  2.0,  together  with  the  contours  of  the  streamwise  velocity  at 
y+  =  7.  The  A  shock  is  consistent  with  the  experimental  Schlieren  image  (Fig.  41). 

In  summary,  the  following  conclusions  can  be  made: 

•  The  LES  does  not  accurately  reproduce  the  mean  surface  pressure.  This  may  be 
attributable  to  the  difference  in  Reynolds  numbers  between  the  LES  and  experiment; 
however,  the  results  of  Kannepalli  et  al  [86]  (Case  3)  at  nearly  the  same  conditions 
as  the  present  case  show  good  agreement  with  the  experimental  surface  pressure  at  a 
comparably  larger  Reynolds  number,  so  the  difference  in  this  case  is  yet  unresolved. 

•  The  LES  predicts  a  mean  separation  length  consistent  with  a  linear  extrapolation  of 
the  experimental  data. 

2.5  Case  5:  Mach  2.3,  a  =  8° 

Gamier  et  al  [101,  102]  simulated  a  Mach  2.3  impinging  shock  at  Res  =  6.0  x  104.  The 
incident  shock  angle  is  32.4°  corresponding  to  a  wedge  angle  a  =  8°.  The  wall  is  adiabatic. 
Experimental  data  of  Deleuze  [103]  and  Laurent  [104]  is  available.  The  flow  conditions  are 
summarized  in  Table  10.  The  overall  inviscid  pressure  ratio  across  the  shocks  P2/P1  =  2.47. 

Details  of  the  computation  are  presented  in  Table  11.  A  structured  grid  LES  algorithm 
was  utilized.  Three  computations  (denoted  by  LES-X  (a)  to  LES-X  (c))  were  performed 
using  the  Mixed  Scale  Model  [105,  106,  107],  and  one  computation  using  MILES  (denoted 
by  LES-M  (d)).  The  inflow  conditions  were  obtained  from  a  separate  computation  of  a  flat 
plate  adiabatic  boundary  layer.  The  inflow  velocity  satisfies  the  compressible  Law  of  the 
Wall  as  indicated  in  Fig.  42  where  HWA  indicates  the  experimental  data  obtained  using  Hot 
Wire  Anemometry. 

The  mean  skin  friction  coefficient  is  shown  in  Fig.  43.  The  computations  using  the 
Mixed  Scale  Model  are  in  close  agreement  with  experiment  and  each  other.  The  MILES 
computation  (Case  D)  also  displays  close  agreement  with  experiment,  and  exhibits  a  smaller 
separation  region  consistent  with  higher  inflow  turbulence  level  near  the  wall  [101,  102], 

The  evolution  of  the  compressible  displacement  thickness  <5i  is  displayed  in  Fig.  44.  All 
computations  show  good  agreement  with  experiment  within  the  interaction  region,  although 
overestimating  5i  in  the  recovery  region  downstream  by  approximately  10%. 

The  mean  streamwise  velocity  profiles  at  x  =  351  mm  (computations)  and  x  =  345  mm 
(experiment)  are  shown  in  Fig.  45  where  LDA  indicates  the  experimental  measurements 
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using  Laser  Doppler  Anemometer.  An  additional  profile  for  Case  A  is  shown  at  x  —  345  mm. 
The  x  location  is  immediately  downstream  of  reattachment  (Fig.  43).  Note  the  velocity  is 
the  ordinate  and  the  distance  from  the  wall  (normalized  by  the  incoming  boundary  layer 
thickness)  is  the  abscissa.  The  offset  between  the  x  locations  for  the  computations  and 
experiment  accounts  for  the  apparent  displacement  of  the  incident  shock  in  the  computations. 
The  overall  agreement  between  computation  and  experiment  is  good. 

The  rms  streamwise  fluctuating  velocity  profiles  at  x  —  326  mm  (computations)  and 
x  =  320  mm  (experiment)  are  shown  in  Fig.  46,  and  at  x  =  351  mm  (computations) 
and  x  =  345  mm  (experiment)  in  Fig.  47.  An  additional  profile  for  Case  A  is  shown  at 
x  =  320  mm  and  x  =  345  mm  in  Figs.  46  and  47,  respectively.  The  profiles  in  Fig.  46 
correspond  to  the  separated  region,  and  the  profiles  in  Fig.  47  are  located  downstream  of 
the  reattachment  of  the  boundary  layer  (Fig.  43).  In  Fig.  46,  the  computed  profiles  display 
general  agreement  with  experiment  (LDA).  In  Fig.  47,  there  is  a  significant  difference  between 
the  HWA  and  LDA  measurements  especially  in  the  middle  of  the  boundary  layer.  The 
computations  show  good  agreement  with  the  LDA  data. 

In  summary,  the  following  conclusions  can  be  made: 

•  The  LES  accurately  predicts  the  displacement  thickness,  skin  friction  coefficient,  and 
velocity  profiles 

•  The  LES  shows  reasonable  agreement  with  the  LDA  measurements  of  rms  streamwise 
turbulent  fluctuations 

•  The  differences  between  the  Mixed  Model  and  MILES  results  are  small 


3  2-D  Expansion-Compression  Corner 

The  expansion-compression  corner  (Fig.  3)  is  a  typical  configuration  to  study  the  response 
of  a  compressible  turbulent  boundary  layer  to  the  effects  of  a  favorable  pressure  gradient, 
convex  streamline  curvature  and  bulk  dilatation.  As  a  supersonic  flow  passes  through  an 
expansion  corner,  the  flow  is  accelerated  through  an  expansion  fan  formed  around  the  corner 
and  the  boundary  layer  becomes  thicker  due  to  the  decrease  in  density.  The  inviscid  flowfield 
immediately  behind  the  expansion  fan  can  be  obtained  from  Prandtl-Meyer  theory.  The 
boundary  layer  experiences  a  relaxation  process  after  the  expansion  fan  and  exhibits  different 
turbulence  features  from  the  incoming  equilibrium  flow.  Several  critical  features  of  this  flow 
configuration  will  be  addressed. 

The  primary  experimental  observation  is  the  increase  of  the  boundary  layer  thickness.  In 
the  measurements  of  Dawson  et  al  [108]  and  Arnette  et  al  [109]  for  a  Mach  3  turbulent  bound¬ 
ary  layer  through  centered  or  gradual  expansions  of  7°  and  14°,  the  time-average  Schlieren 
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images  indicate  that  the  boundary  layer  thickness  increases  approximately  150%  across  the 
7°  expansion  and  200%  across  the  14°  expansion  (Fig.  48).  It  has  been  demonstrated  by 
Morkovin  [110]  that  the  induced  bulk  dilatation  serves  to  increase  the  boundary  layer  thick¬ 
ness.  The  instantaneous  Filtered  Rayleigh  Scattering  (FRS)  visualization  (Fig.  49)  shows 
that  the  structures  undergo  an  increase  in  both  scale  and  angular  orientation  behind  the 
expansion.  Arnette  et  al  [109]  suspected  that  the  increase  in  structure  angle  is  a  kinematic 
effect  associated  with  the  acceleration  of  the  bottom  of  the  structure  through  the  inclined 
expansion  region  before  the  top  of  the  structure. 

The  reduction  of  turbulence  is  another  significant  feature  of  this  flow  configuration.  This 
flow  configuration  experiences  the  suppression  of  turbulence  due  to  the  combined  effects  of 
the  favorable  pressure  gradient,  the  convex  streamline  curvature  and  the  bulk  dilatation.  The 
reduction  in  the  turbulence  may  be  strong  enough  to  relaminarize  the  flow.  Early  work  by 
Morkovin  [110]  found  that  longitudinal  turbulence  intensities  and  wall  shear  stress  decreased 
significantly  due  to  a  favorable  pressure  gradient  and  a  bulk  dilatation. 

Thomann  [111]  isolated  the  effect  of  streamline  curvature  on  the  heat  transfer  rate  at 
the  wall  for  a  Mach  2.5  and  20°  of  convex  curvature  turbulent  boundary  layer  by  placing 
appropriately  shaped  bodies  in  the  free  stream  above  convex  and  concave  surfaces  to  elim¬ 
inate  pressure  gradients.  The  heat  transfer  rate  at  the  wall  decreased  approximately  20%, 
indicating  decreased  turbulent  mixing  between  the  boundary  layer  and  the  freestream. 

Lewis  et  al  [112]  studied  a  Mach  4  turbulent  boundary  layer  subject  to  isolated  adverse 
and  favorable  pressure  gradients.  The  experiments  showed  that  under  a  favorable  pressure 
gradient  the  ability  of  the  turbulence  to  transport  momentum  to  the  wall  was  markedly 
decreased,  which  was  largely  due  to  the  effect  of  the  strong  expansion  on  intermittency. 
In  the  analytical  analysis,  a  compressible  pressure  gradient  parameter  (3k  =  5*K / rw(dp / dx) 
(where  5*K  is  the  kinematic  displacement  thickness)  suggested  by  Alber  and  Coats  [113]  was 
used  as  the  dominant  parameter.  The  boundary  layer,  when  subjected  to  a  large  negative 
(3k,  was  shown  to  be  poorly  characterized  by  the  Law  of  the  Wake  and  was  found  to  relax  to 
an  equilibrium-like  profile  after  fifty  initial  boundary  layer  thickness.  A  similar  parameter 
suggested  by  Narasimha  and  Sreenivasan  [114]  will  be  discussed  later  on  as  an  important 
criterion  for  the  flow  relaminarization. 

So  and  Mellor  [115,  116]  performed  a  series  of  experiments  providing  detailed  turbulence 
measurements  in  boundary  layers  perturbed  by  convex  curvature.  The  parameter  8/R  was 
kept  constant  at  approximately  0.08  and  0.12,  respectively,  over  a  distance  of  approximately 
505o  where  5  and  <5o  are  local  and  incoming  boundary  layer  thickness,  respectively,  and  R  is 
the  radius  of  the  curvature.  These  experiments  showed  that  the  shear  stress  for  y/5  >  0.5 
collapsed  to  virtually  zero  within  a  few  boundary-layer  thicknesses.  The  absence  of  shear 
stress  in  the  outer  layer  means  there  is  no  turbulence  production  mechanism  and  consequently 
the  turbulence  there  decays  as  it  moves  downstream.  The  most  interesting  aspect  was  a 


RTO-TR-AVT-007-V3 


11  -23 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


ORGANIZATION 


change  of  scale.  After  a  sustained  region  of  convex  curvature,  —u'v'/u/  scaled  with  y/R 
rather  than  y/5. 

Gillis  et  al  [117,  118]  studied  the  effect  of  a  convex  curvature  on  turbulent  boundary  layers. 
Two  experiments  were  performed,  with  5q/ R  approximately  0.1  and  0.05,  respectively.  The 
curvature  persisted  for  about  20  J0  and  40  J0,  respectively  and  the  total  turning  angle  in  both 
cases  was  90°.  Thereafter,  the  curvature  was  suddenly  removed  and  the  boundary  layer  was 
allowed  to  relax  on  a  flat  plate  with  zero  pressure  gradient.  The  experiment  found  that  the 
flow  showed  a  very  slow  recovery;  the  skin  friction  coefficient  showed  no  sign  of  recovery  for 
about  40ho  downstream,  although  this  recovery  behavior  appears  to  scale  with  x/R  rather 
than  x/5q.  This  slow  response  is  a  direct  result  of  the  almost  total  absence  of  shear  stress  in 
the  outer  layer,  similar  to  the  observations  by  So  and  Mellor  [115,  116].  Those  observations 
led  Gillis  and  Johnston  to  propose  a  two-layer  model  for  the  perturbed  boundary  layer, 
consisting  of  an  inner  active  layer  where  the  turbulence  is  still  in  production  and  an  outer 
layer  containing  nearly  isotropic  “debris”  from  the  thick,  turbulent  boundary  layer  upstream 
of  curvature.  The  recovery  is  dominated  by  the  slow  growth  of  the  inner  active  layer  and 
the  newly  created  turbulence  in  the  inner  layer  retains  little  memory  of  the  upstream  history 
[119]. 

Dussauge  and  Gaviglio  [120]  studied  the  behavior  of  a  boundary  layer  as  it  passes  through 
an  expansion  produced  by  a  sudden  deflection  of  the  wall  at  a  sharp  edge.  The  flow  expe¬ 
rienced  a  short  region  of  favorable  pressure  gradient,  as  well  as  the  stabilizing  effects  of 
convex  curvature  and  bulk  dilatation.  Dussauge  and  Gaviglio  measured  the  longitudinal 
turbulence  intensity  and  the  static  temperature  fluctuation  intensity  and  they  found  that 
u'2  was  reduced  to  half  its  upstream  value.  However,  Morkovin’s  “strong  Reynolds  analogy” 
[121]  is  still  valid  supported  by  the  fact  that  the  temperature- velocity  correlation  coefficient 
remained  constant  at  a  value  of  -0.8.  The  dilatation  appeared  to  be  responsible  for  about 
two-thirds  of  the  observed  reduction  in  u'2,  with  the  pressure-strain  term  contributing  to 
the  remaining  third.  The  choice  for  the  pressure-strain  term  model  was  not  crucial;  three 
different  models  all  gave  about  the  same  result. 

Dussauge  and  Gaviglio  [122]  investigated  a  12°  centered  expansion  of  a  Mach  1.76  equi¬ 
librium  turbulent  boundary  layer  at  Reg  =  5000  and  Jo  —  10  mm  where  Jo  is  incoming 
boundary  layer  thickness.  Experiments  showed  that  the  mean  velocity  profiles  downstream 
of  the  expansion  initially  displayed  a  thick  sublayer  region  possessing  a  larger  normal  gradi¬ 
ent  of  mean  streamwise  velocity  than  the  equilibrium  profiles,  with  no  apparent  logarithmic 
region  until  approximately  10Jo  downstream  of  the  corner.  The  reduction  of  longitudinal 
turbulent  intensities  across  the  entire  boundary  layer  was  observed  in  the  experiments.  The 
reductions  were  substantial  closer  to  the  wall  while  they  became  smaller  for  y/5  >  0.7  (where 
J  is  the  local  boundary  layer  thickness).  Downstream  measurements  showed  the  near- wall 
region  to  re-establish  turbulence  intensities  well  below  the  incoming  levels  more  quickly  (10Jq 
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downstream  of  the  corner)  than  the  outer  part  of  the  boundary  layer,  where  the  evolution 
was  slow.  This  rapid  recovery  of  the  region  near  the  wall  led  Dussauge  and  Gaviglio  to  pro¬ 
pose  that  a  new  internal  layer  formed  downstream  of  the  expansion  and  that  the  incoming 
boundary  layer  had  been  relaminarized. 

Rapid  Distortion  Theory  was  used  to  isolate  the  effects  of  bulk  dilatation  by  neglecting  all 
the  triple-product  terms  and  the  dissipation  terms.  The  calculation  reproduced  the  reduction 
in  longitudinal  turbulence  intensity  very  well  in  the  outer  portions  of  the  boundary  layer, 
but  provided  a  poor  prediction  near  the  wall.  This  indicated  that  the  bulk  dilatation  played 
a  relatively  less  important  role  near  the  wall.  It  also  indicated  that  while  Rapid  Distortion 
Theory  may  be  capable  of  describing  the  initial  response  to  a  sudden  perturbation,  it  cannot 
in  principle  describe  the  whole  relaxation  process. 

Smith  and  Smits  [123]  studied  a  20°  centered  expansion  of  a  Mach  2.84  boundary  layer 
at  Reg  =  77600  and  5q  =  26  mm.  Measurements  at  3.5<5o  downstream  of  the  expansion 
corner  compared  favorably  with  those  of  Dussauge  and  Gaviglio  [122]  in  which  no  evident 
logarithmic  region  in  the  mean  velocity  profile  and  both  the  reduced  longitudinal  turbulence 
intensity  and  Reynolds  stress  were  observed.  A  prominent  observation  was  that  the  mass 
flux  fluctuation  profile  remained  unchanged  through  the  expansion.  This  raised  some  doubts 
about  the  relaminarization. 

The  relaminarization  of  a  turbulent  flow  when  it  passes  an  expansion  corner  has  been 
reported  by  many  investigators  [110,  122,  123,  124,  125,  126,  127,  128].  In  particular, 
Narasimha  and  Sreenivasan  [114]  proposed  that  relaminarization  occurs  when  the  ratio  of 
the  pressure  gradient  to  a  characteristic  Reynolds  stress  gradient  ( 5dp/dx/rw )  is  larger  than 
75,  where  5  is  the  local  boundary  layer  thickness  and  rw  is  the  wall  shear  stress  in  the  bound¬ 
ary  layer  just  upstream  of  the  expansion  corner.  Taking  dp/dx  ~  A  p/6,  this  correlation  is 
simplified  to  Ap/ro  >  75  [129],  where  A p  is  the  pressure  drop  across  the  expansion  fan.  In 
the  study  by  Smith  and  Smits  [123],  although  this  published  criterion  for  relaminarization 
had  been  met,  caution  was  urged  in  using  the  term. 

Arnette  et  al  [109,  125,  130]  used  Filtered  Rayleigh  Scattering  to  investigate  the  effects 
of  expansions  on  the  Mach  3  turbulent  boundary  layer  subjected  to  four  expansion  regions 
(centered  and  gradual  expansions  of  7°  and  14°).  An  examination  of  the  compressible  vortic- 
ity  transport  equation  and  estimates  of  the  perturbation  impulses  attributable  to  streamline 
curvature,  acceleration  and  dilatation  both  confirm  dilatation  to  be  the  primary  stabilizer. 

Experiments  showed  that  the  large-scale  structures  of  the  outer  portions  of  the  boundary 
layer  maintained  their  identity  and  increased  in  scale  across  the  expansions  and  the  structure 
angle  also  appeared  to  increase,  suggesting  small-scale  motions  were  quenched  and  conse¬ 
quently  resulting  in  a  more  intermittent  boundary  layer,  consistent  with  the  observations  by 
Lewis  et  al  [112],  However,  the  large-scale  structures  weakened  downstream  of  the  expansion 
regions  supported  by  the  reduction  of  the  Reynolds  shear  stress  which  is  mainly  associated 
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with  large-scale  structures. 

Turbulence  levels  were  observed  to  decrease  across  the  expansion  and  the  reductions 
increased  in  severity  as  the  wall  is  approached  since  the  turbulence  activity  is  associated 
with  small-scale  motions  near  the  wall  which  are  significantly  quenched.  This  observation 
seems  to  agree  with  the  results  of  Dussauge  and  Gaviglio  [122]  and  Smith  and  Smits  [123], 
where  sharp  reductions  in  near  wall  turbulence  were  encountered  across  the  expansions. 
Additionally,  the  severity  of  the  reductions  increased  with  expansion  angle.  Downstream  of 
the  two  14°  expansions,  the  reductions  are  more  severe  and  reverse  transition  is  indicated 
by  sharp  reductions  in  turbulent  kinetic  energy  levels  and  a  change  in  sign  of  the  Reynolds 
shear  stress,  while  the  two  7°  expansions  are  not  strong  enough  to  cause  complete  reversion. 

For  all  of  the  expansions,  reductions  in  the  Reynolds  shear  stress  levels  occur  more 
rapidly  than  decreases  in  the  turbulence  levels.  The  Reynolds  shear  stress  in  part  of  the  top 
half  of  the  boundary  actually  changes  sign,  indicating  that  the  usual  extraction  of  energy 
from  the  mean  flow  by  the  turbulence  has  ceased  and  the  turbulent  motions  are  decaying. 
Furthermore,  reductions  in  the  turbulent  transport  of  turbulent  kinetic  energy  in  the  normal 
direction  occur  more  rapidly  than  reductions  for  the  turbulent  transport  in  the  streamwise 
direction,  as  indicated  by  the  fact  that  streamwise  skewness  and  flatness  profiles  were  not 
altered  significantly,  but  the  normal  skewness  was  reduced  and  the  normal  flatness  was 
decreased. 

The  recovery  process  is  not  only  associated  with  the  regions  of  the  boundary  layer  [122], 
but  also  with  the  quantities  measured.  Experiment  showed  that  the  mean  velocity  profiles 
would  recover  the  equilibrium  profile  more  quickly  than  the  turbulence  profiles.  Indications 
of  recovery  appear  nearer  the  expansion  regions  for  the  7°  expansions  than  for  the  14° 
expansions.  Bradshaw  [131]  suggested  that  a  suitable  response  time  for  the  stress-containing 
eddies  is  the  ratio  of  the  turbulence  energy  k  to  its  rate  of  production  (f  =  k/(u'v'  |^),  where 

the  stress  ratio  a  =  and  the  typical  value  is  0.3).  A  corresponding  relaxation  length  X 
is  u  times  this  response  time.  In  the  middle  of  the  layer,  using  a  1  /7th  power  law,  u/^  ~  3 
and  therefore  X  &  105o-  It  is  generally  true  that  the  outer  layer  responds  very  slowly  to 
a  disturbance.  As  the  wall  is  approached,  X  tends  to  zero;  that  is,  the  flow  adjusts  very 
quickly  in  the  inner  region.  The  relaxation  length  also  depends  on  the  quantity  measured. 
This  occurs  largely  because  of  the  increasing  sensitivity  of  the  turbulence  terms  as  they 
increase  in  order.  For  example,  the  shear  stress  may  take  longer  than  the  mean  velocity  to 
relax,  partly  because  it  becomes  self-preserving  at  roughly  the  same  rate  as  du/dy. 

Dawson  et  al  [108]  studied  the  response  of  a  Mach  3  fully  developed  compressible  turbu¬ 
lent  boundary  layer  at  Reg  =  2.5  x  104  to  centered  and  gradual  (R/50  =  50)  expansions,  both 
of  7°  and  14°  deflection.  High-frequency-response  miniature  pressure  transducers  were  used 
to  acquire  multi-point  fluctuating  surface  pressure  measurements  in  the  same  equilibrium 
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and  perturbed  boundary  layers.  Normalized  power  spectra  showed  the  pressure  fluctuations 
to  be  much  more  concentrated  at  low  frequencies  just  downstream  of  the  expansions  rel¬ 
ative  to  the  pre-expansion  flow.  The  elevated  low-frequency  levels  were  accompanied  by 
sharp  decreases  at  high  frequencies.  This  indicated  that  small-scale  motions  were  quenched 
essentially  immediately  after  the  expansion. 

Lu  and  Chung  [132]  proposed  a  scaling  law  to  identify  a  downstream  influence  of  the 
expansion  corner  based  on  the  mean  surface  pressure  distribution.  The  downstream  influence 
of  the  corner  is  estimated  as  the  horizontal  distance  from  the  corner  to  the  intersection  of  the 
tangent  through  the  downstream  pressure  data  with  the  inviscid  downstream  pressure.  The 
scaling  parameter  is  chosen  to  be  hypersonic  similarity  parameter  K  =  where  M, x 

is  incoming  Mach  number  and  a  is  the  corner  deflection  angle  in  radians.  The  downstream 
influence  distance  increases  with  K.  This  scaling  law  is  validated  for  Mach  1.76  to  Mach  8 
for  K  up  to  1. 

Chung  and  Lu  [133]  study  experimentally  the  effect  of  the  expansion  on  surface  pressure 
fluctuations.  Surface  pressure  fluctuations  of  Mach  8  turbulent  flow  past  a  2.5°  and  a  4.25° 
expansion  corner  maintained  a  Gaussian  distribution  but  were  severely  attenuated  by  the 
expansion  process.  The  surface  pressure  fluctuations  four  boundary-layer  thickness  down¬ 
stream  of  the  corner  damped  to  only  37%  and  21%  of  the  incoming  value  for  the  2.5°  and 
a  4.25°  corners.  The  pressure  fluctuations  did  not  recover  to  the  level  corresponding  to  an 
equilibrium  turbulent  flow  even  though  the  mean  pressure  downstream  reached  the  inviscid 
value  within  four  to  six  boundary  layer  thicknesses.  The  space-time  correlation  study  of  the 
surface  pressure  fluctuations  showed  that  the  fluctuations  were  convected  with  a  velocity 
comparable  to  that  on  a  flat  plate  and  these  fluctuations  maintained  their  identities  longer 
for  stronger  expansions. 

Chung  [134]  performed  experiments  for  the  incoming  flow  past  5°,  10°  and  15°  expansion 
corners  at  M  =  1.28  and  the  unit  Reynolds  number  of  2.57  x  107m-1.  The  study  showed 
that  the  downstream  influence  for  a  lower  supersonic  expansion-corner  flow  at  Mm  =  1.28 
is  insensitive  to  the  corner  deflection  angle  (a  =  5°,  10°  and  15°).  The  downstream  influ¬ 
ence  was  measured  by  the  mean  surface  pressure  and  the  surface  pressure  fluctuation.  The 
downstream  influence  is  less  than  one  boundary-layer  thickness  based  on  the  mean  sur¬ 
face  pressure  distribution.  This  agrees  with  the  discussion  of  Narasimha  and  Sreenivasan 
[135].  The  surface  pressure  fluctuations  are  normalized  by  upstream  dynamic  pressure.  The 
distributions  indicated  that  the  intensity  of  pressure  fluctuations  decreases  and  reaches  the 
minimum  downstream  of  expansion  corners.  The  final  equilibrium  pressure  fluctuation  levels 
are  obtained  with  three-  to  five-boundary  layer  thickness,  which  indicates  a  larger  interac¬ 
tion  region  than  obtained  from  the  mean  surface  pressure  distribution.  Thus  the  interaction 
region  scaling  of  the  expansive  flows  may  be  further  refined  based  on  the  distribution  of 
surface  pressure  fluctuation. 
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Zheltovodov  et  al  [24]  performed  a  combined  experimental  and  numerical  study  of  a 
supersonic  turbulent  flow  over  an  expansion  corner  followed  by  a  compression  corner  at 
different  Mach  numbers  and  deflection  angles.  The  effects  of  the  expansion  corner  on  the 
separation  properties  were  thoroughly  studied.  It  was  shown  that  the  basic  characteristics 
of  the  flow  regimes  are  1)  attached  flow  for  small  values  of  a,  2)  the  appearance  of  a  local 
separation  zone  in  the  compression  corner  with  a  free  separation  point  at  moderate  values  of 
a,  3)  formation  of  large-scale  separated  flow  with  a  fixed  separation  point  when  a  becomes 
large  enough.  A  range  of  the  dimensionless  separation  length  is  plotted  vs  Reynolds  number 
in  Fig.  50  where  the  characteristic  length  (Lc)  is  defined  in  (1)  with  S  and  evaluated 
downstream  of  the  expansion  fan  and  upstream  of  the  compression  corner.  As  shown  by 
Zheltovodov  et  al  [23,  24],  the  increase  of  the  velocity  profile  ” fullness”  in  the  boundary 
layer  after  the  expansion  fan  can  increase  its  resistance  to  separation  compared  to  an  isolated 
compression  corner  case  (see  Fig.  7).  The  monotonic  increase  in  the  separation  length  with 
the  decrease  in  Reynolds  number  is  observed  in  Figs.  7  and  50.  This  can  be  explained  by  a 
decrease  in  the  effective  viscosity  near  the  wall  due  to  the  fact  that  the  turbulent  fluctuations 
of  the  expansion  waves  become  more  suppressed  as  the  Reynolds  number  decreases.  This 
is  in  accordance  with  the  conclusion  of  Goldfeld  et  al  [128]  regarding  the  Reynolds  number 
effects  on  relaminarization  of  the  turbulent  boundary  layer  behind  a  convex  corner. 

Zheltovodov  et  al  [136]  analyzed  the  relaxation  processes  in  the  disturbed  turbulent 
boundary  layer  in  the  vicinity  of  expansion-compression  corner  configuration  at  different 
surface  inclination  angles  8,  25  and  45  degrees  on  the  basis  of  the  velocity  profiles  and 
surface  skin  friction  measurements  at  Mach  3  and  4.  The  appearance  of  the  distinct  Gortler 
type  vortices  has  been  ascertained  by  the  surface  pattern  visualization  in  the  vicinity  of  the 
45  deg  configuration.  A  comparison  of  these  data  with  similar  results  for  the  compression- 
expansion  corner  configuration  at  the  same  surface  inclination  angles  has  been  performed 
to  determine  the  complex  nature  of  the  dependence  of  the  relaxation  process  on  the  shock 
wave  and  expansion  fan  positions. 

Zheltovodov  et  al  [37,  39,  40,  41]  performed  measurements  of  surface  heat  transfer  dis¬ 
tribution  in  the  vicinity  of  25  degree  expansion-compression  corner  Mach  2,  3  and  4.  As 
shown  by  the  experiments,  the  decrease  in  Mach  number  caused  an  increase  in  the  length  of 
the  separation  zone  in  the  compression  corner  after  the  expansion  fan  and  the  heat  transfer 
level  in  it,  as  well  as  an  increase  of  the  heat  transfer  level  on  the  plate  downstream  of  the 
attachment  line.  The  heat  transfer  distribution  for  such  configuration  has  been  compared 
with  similar  data  for  the  25  deg  compression-expansion  corner  case.  As  shown  on  the  basis 
of  analysis  of  the  hot-wire  measurements  of  the  mass  flow,  density  and  velocity  fluctuations, 
the  higher  heat  transfer  level  downstream  of  the  compression/expansion  corner  comparing 
with  the  expansion/compression  configuration  is  associated  with  a  higher  turbulence  level 
downstream  of  the  disturbances  in  the  first  case.  Zheltovodov  et  al  [39,  40]  and  Borisov 
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et  al  [48]  demonstrated  the  possibility  of  prediction  of  the  heat  transfer  distribution  in  the 
vicinity  of  such  configurations  on  a  basis  of  the  asymptotic  turbulent  boundary  layer  the¬ 
ory  by  Kutateladze-Leontiev  in  the  framework  of  the  ”  new-layer”  concept  with  the  use  of 
generalized  experimental  turbulence  data. 

Zhcltovodov  et  al  [24],  Borisov  et  al  [48,  137],  Horstman  and  Zhcltovodov  [64,  138] 
used  the  experimental  data  for  the  8,  25  and  45  deg  expansion-compression  corner  flows 
at  Mach  3  and  4  to  validate  RANS  simulations  with  different  turbulence  models  (k  —  e  by 
Jones  and  Launder  and  its  two  layer  modification  by  Rodi,  q  —  oo  by  Coakley).  Additional 
review  of  these  studies  has  been  performed  by  Zheltovodov  [3,  47].  It  was  concluded  that  the 
mean  flowheld  structure,  surface  pressure  and,  sometimes,  skin  friction  distribution  can  be 
accurately  predicted.  Nevertheless,  surface  heat  transfer  distribution  cannot  be  accurately 
predicted  by  the  the  turbulence  models  used. 

In  the  following  section,  an  LES  of  the  expansion-compression  corner  is  presented.  Mean 
quantities  represent  combined  time  and  spanwise  averages. 

3.1  Case  1:  Mach  3,  a  =  25° 

Knight  et  al  [139]  simulated  a  Mach  2.88  expansion  compression  corner  at  Res  =  2  x  104 
and  adiabatic  wall  conditions.  Experimental  data  of  Zheltovodov  et  al  [23,  24,  37,  39,  40, 
41,  136,  138]  is  available  for  a  =  25°  at  =  2.9  and  Res  =  4.1  x  104  to  1.95  x  105.  The 
flow  conditions  are  shown  in  Table  12.  The  inviscid  pressure  ratio  across  the  expansion  is 
P2/P1  =  0.098,  and  the  Mach  number  downstream  of  the  expansion  is  M2  =  4.59.  Across 
the  subsequent  compression,  the  pressure  ratio  is  Ps/p-2  =  8.60  and  the  downstream  Mach 
number  M3  =  2.43. 

Details  of  the  computation  are  presented  in  Table  13.  An  unstructured  grid  LES  algo¬ 
rithm  using  the  MILES  model  [88]  was  utilized.  The  inflow  profile  is  the  same  as  Case  1  of 
the  2-D  compression  corner. 

The  computed  flowheld  structure  is  shown  in  Figs.  51  and  52  which  display  the  mean 
static  pressure  and  streamlines  (S  is  separation,  A  is  attachment).  The  how  expands  around 
the  hrst  corner,  and  recompresses  at  the  second  corner.  The  shock  wave  turbulent  boundary 
layer  at  the  second  corner  is  sufficiently  strong  to  separate  the  boundary  layer  as  evident  in 
Fig.  52.  The  howheld  structure  is  in  good  agreement  with  the  results  of  Zheltovodov  and 
his  colleagues  [3,  23,  24,  37,  40,  41,  136]  which  are  shown  qualitatively  in  Fig.  3. 

Mean  velocity  prohles  in  the  x-direction  are  shown  in  Fig.  53  at  x  =  26  and  x  =  6 6,  where 
x  is  measured  from  the  inflow  along  the  direction  of  the  inflow  freestream  velocity  (Fig.  51). 
The  abscissa  is  the  component  of  velocity  locally  parallel  to  the  wall,  and  the  ordinate  is  the 
distance  measured  normal  to  the  wall.  The  hrst  profile  is  upstream  of  the  expansion  corner 
(which  is  located  at  x  —  4 6).  The  second  profile  is  downstream  of  the  expansion  fan  and 
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upstream  of  the  separation  point.  The  computed  mean  velocity  profile  at  the  first  location 
is  slightly  fuller  than  the  experiment.  This  is  consistent  with  the  experimentally  observed 
dependence  of  the  exponent  n  in  the  power-law  U /[Too  =  (y/^)1^  011  the  Reynolds  number. 
The  second  profile  shows  a  significant  acceleration  of  the  flow  in  the  outer  portion  of  the 
boundary  layer  due  to  the  expansion. 

The  surface  pressure  is  displayed  in  Fig.  54.  The  surface  pressure  drops  rapidly  at  the 
first  corner.  The  shock  boundary  layer  interaction  at  the  second  corner  causes  separation  of 
the  boundary  layer  and  a  concomitant  plateau  in  the  surface  pressure  [67].  The  experiments 
display  an  increase  in  the  size  of  the  pressure  plateau  region  with  decreasing  Reynolds  num¬ 
ber.  The  computed  results  for  Res  =  2  x  104  are  in  good  agreement  with  the  experimental 
data  at  the  lowest  Reynolds  number  (Res  =  4.1  x  104)  for  the  location,  extent  and  magnitude 
of  the  pressure  plateau.  Also,  the  shape  of  the  experimental  pressure  plateau  shows  little 
variation  for  Res  <  6.8  x  104,  thus  suggesting  that  the  computed  pressure  plateau  region 
(for  Res  =  2  x  104)  is  accurate.  The  computed  recovery  of  the  surface  pressure  is  more  rapid 
than  in  the  experiment,  however. 

The  computed  and  experimental  mean  skin  friction  coefficient  Cf  =  Tw/^p^U^  are  shown 
in  Fig.  55.  The  skin  friction  coefficient  increases  rapidly  downstream  of  attachment.  The 
computed  results  at  Res  =  2  x  104  are  in  close  agreement  with  the  experimental  data  at 
Res  =  8.0  x  104  and  1.94  x  105  in  the  region  downstream  of  reattachment. 

Zhcltovodov  and  Schiilein  [23]  have  correlated  experimental  data  for  the  scaled  sepa¬ 
ration  length  Lsep/Lc  defined  by  (2)  for  the  expansion-compression  corner  where  6  is  the 
incoming  boundary  layer  thickness  (upstream  of  the  expansion  corner),  P2  is  the  pressure 
after  the  shock  in  inviscid  flow,  and  ppi  is  the  plateau  pressure  from  the  empirical  formula 
of  Zukoski  [25]  ppi  =  pe(^Me  +  1)  where  pe  and  Me  are  the  static  pressure  and  freestream 
Mach  number  upstream  of  the  compression  corner  and  downstream  of  the  expansion  fan.  I11 
the  computation,  the  location  is  taken  to  be  x  =  6<5.  The  values  of  Me  and  P2  have  been 
computed  using  inviscid  theory.  Also,  Rese  =  1-8  x  104  for  LES  (Rese  =  PeUeSe/pe,  where 
pe,Ue  and  pe  are  computed  using  inviscid  theory).  The  experimental  data  correlation  and 
the  computed  result  [139]  for  the  scaled  separation  length  is  shown  in  Fig.  50.  The  computed 
value  is  consistent  with  a  linear  extrapolation  of  the  experimental  data. 

I11  summary,  the  following  conclusions  can  be  made: 

•  The  LES  accurately  reproduces  the  mean  surface  pressure  and  rapid  recovery  of  the 
skin  friction  downstream  of  attachment  to  a  level  close  to  the  undisturbed  upstream 
in  agreement  with  experiment 

•  The  LES  predicts  a  mean  separation  length  consistent  with  a  linear  extrapolation  of 
the  experimental  data 
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4  3-D  Single  Fin 

The  3-D  single  fin/plate  configuration  (Fig.  4)  is  a  standard  CFD  validation  case  for  a 
swept  shock  wave  turbulent  boundary  layer  interaction.  An  equilibrium  supersonic  turbulent 
boundary  layer,  developed  along  a  plate,  approaches  a  semi-infinite  fin  of  angle  a ,  attached 
normally  to  the  plate.  The  deflection  of  the  flow  by  the  fin  generates  a  swept  shock  system. 
The  increase  in  pressure  through  the  shock  is  diffused  on  the  plate  and  a  disturbed  flow 
pattern  is  observed  for  a  considerable  distance  both  upstream  and  downstream  of  the  shock 
position  predicted  for  inviscid  flow.  If  the  shock  is  strong  enough  to  cause  the  boundary  layer 
to  separate,  the  topology  of  the  flow  changes  significantly.  The  intensity  of  the  interaction 
depends  on  the  flow  conditions  (Mach  number,  fin  angle),  and  on  the  geometry  of  the  fin 
(swept  or  not). 

In  this  section  we  review  the  most  critical  features  of  the  single  fin  flow  field.  A  more 
detailed  review  has  been  done  by  Panaras  [140].  Additional  reviews  are  presented  in  Settles 
and  Dolling  [141,  142],  Zheltovodov  and  Schiilein  [143],  Bogdonoff  [144],  Zheltovodov  [3] 
and  Knight  [1].  In  general  the  interaction  domain  in  these  types  of  flows  is  quasi-conical, 
i.e.,  it  grows  almost  linearly  in  the  downstream  direction.  Various  manifestations  of  the 
conical  nature  of  the  flow  have  been  observed  experimentally,  even  in  the  early  70’s,  when  the 
available  experimental  techniques  were  limited.  Nevertheless,  it  was  very  recently  that  strong 
evidence  was  provided  for  the  existence  of  quasi-conical  vortices  in  swept  shock/boundary 
layer  interactions. 

The  existence  of  pressure  gradients  in  the  transverse  direction  as  well  as  along  the  external 
streamline  direction  is  the  essential  characteristic,  which  distinguishes  a  3-D  flow  from  a  2- 
D  or  axisymmetric  one.  The  boundary  layer  responds  to  the  transverse  gradient  by  the 
development  of  a  crossflow  or  secondary  flow.  Though  the  crossflow  may  be  small  near  the 
outer  edge  of  the  boundary  layer,  it  may  be  substantially  close  to  the  surface  where  the 
momentum  deficit  is  large.  Thus  the  velocity  vectors  close  to  the  surface  can  progressively 
rotate  resulting  in  a  “skewing”  of  the  boundary  layer.  When  the  distance  from  the  surface 
tends  to  zero  the  velocity  vectors  reach  a  limiting  direction,  which  is  co- linear  with  the  skin 
friction  vectors.  At  the  same  time  the  streamlines  tend  to  a  limiting  position,  which  is  also 
a  trajectory  of  the  skin  friction  lines.  For  this  reason  they  are  called  “limiting  streamlines” . 
In  general,  the  limiting  streamlines  in  3-D  flows  follow  paths  that  are  different  in  direction 
from  the  external  streamlines.  Limiting  streamlines  originate  at  nodal  points  of  attachment, 
and  after  circumscribing  the  body  surface  disappear  into  nodal  points  of  separation. 

In  the  case  of  the  fin/plate  configuration  the  footprint  of  the  flow  (turbulent  or  laminar)  on 
the  surface  of  the  plate  consists  also  of  a  separation  line,  lying  ahead  of  the  inviscid  position 
of  the  shock,  and  of  a  reattachment  line  close  to  the  corner.  Token  [145]  proposed  the  first 
flow  model  of  the  fin-plate  configuration  in  1974.  For  explaining  the  high  heat-transfer  peak 
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measured  on  the  flat  plate  near  the  root  of  the  fin,  he  suggested  that  a  separation  conical 
vortex  appears  between  the  separation  and  the  reattachment  lines.  Kubota  and  Stollery 
[146]  improved  the  flow  model  of  Token  based  on  their  own  experiments.  The  separated  flow 
in  their  model  is  characterized  by  two  counter-rotating  vortices  -  a  tight,  vigorous,  roughly 
circular  one  in  the  corner  with  a  weak,  very  elongated  one  above  it  (Fig.  56).  Kubota 
and  Stollery,  in  addition  to  the  separation  vortices,  investigated  the  structure  of  the  shock 
system.  Of  paramount  importance  for  the  development  of  the  flow  model  is  the  appearance 
or  absence  of  the  A-foot  bifurcation  of  the  fin  shock  at  the  region  of  its  interaction  with  the 
boundary  layer,  exactly  as  it  happens  in  a  strong  2-D  shock- wave /boundary  layer  interaction. 
For  resolving  this  question  Kubota  and  Stollery  applied  the  vapour-screen  technique  to  their 
Mach  3  tests.  According  to  their  vapour-screen  pictures,  when  the  angle  of  the  fin  is  smaller 
than  the  one  required  for  the  appearance  of  separation,  there  exists  no  A-shock,  but  for  larger 
angles,  when  the  corresponding  oil-flow  picture  suggests  a  separated  flow,  there  is  evidence 
of  the  shock  splitting  into  a  A-shape  near  the  edge  of  the  boundary  layer. 

In  the  early  1990’s  non-intrusive  experimental  techniques  were  available  that  produce  very 
clear  pictures  of  the  cross-section  of  the  flowfield  in  the  normal  to  the  shock  direction,  and 
provide  quantitative  data,  such  as  the  skin  friction  distribution,  which  are  suitable  for  direct 
comparison  with  computational  predictions.  Also,  the  efficiency  of  computers  permitted 
the  use  of  high-resolution  grids  and  more  advanced  codes.  The  new  numerical  simulations 
verified  all  the  critical  elements  of  the  flow  structure  of  a  swept  shock/turbulent-boundary- 
layer  interaction,  some  of  which  up  to  then  were  objects  of  speculation. 

An  important  consequence  of  the  quasi-conical  nature  of  the  fin/plate  flows  is  the  pre¬ 
sumption  that  their  features  can  be  projected  upon  the  surface  of  a  sphere  whose  origin  is 
the  conical  origin.  Zubin  and  Ostapenko  [147]  and  Alvi  and  Settles  [148]  have  demonstrated 
this  feature  of  the  conical  flow  by  using  conical  shadowgraphy.  Focusing  a  light  beam  at  the 
origin  of  the  approximately  conical  flowfield  and  aiming  it  such  that  the  resulting  conical 
light  beam  coincided  with  the  rays  of  the  swept  interaction,  they  obtained  clear  pictures 
of  the  flowfield  of  fin/plate  interactions.  Alvi  and  Settles  [149]  combined  their  conical  flow 
images  with  previous  wall  pressure  and  skin-friction  measurements  to  construct  a  physical 
flowfield  model.  An  example  is  shown  in  Fig.  57  for  a  M, x  =  3,  a  =  16°  test  case.  The  cross- 
section  of  the  flow  is  normal  to  the  shock.  It  is  seen  that  on  top  of  the  separation  bubble 
lies  a  well-bifurcated  shock.  From  the  shock  triple-point  a  shear  layer  emanates  and  moves 
towards  the  corner.  Below  the  flowfield  map  the  corresponding  surface  distributions  of  static 
pressure  and  skin  friction  are  plotted  in  the  same  conical-angle  scale.  Alvi  and  Settles  [149] 
observed  that  the  inviscid  air  processed  between  the  triple-point  and  the  separation  vortex 
is  curved  downwards  and  impinges  upon  the  flat  plate.  This  is  accomplished  by  reflected 
Prandtl-Meyer  expansion  and  compression  fans.  The  authors  call  it  an  “impinging-jet”  and 
attribute  to  it  the  peak  heating,  pressure  and  skin  friction,  observed  in  the  corner  region. 
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Clear  3-D  images  of  the  swept  shock/turbulent-boundary-layer  interaction  were  provided 
by  Panaras  [150],  who  calculated  a  =  3.0,  a  =  20°  fin/plate  flow,  using  a  RANS 
code  based  on  the  Baldwin  and  Lomax  [151]  algebraic  turbulence  model.  Fig.  58  from 
Panaras  [150]  includes  all  the  critical  elements  of  the  swept  shock/turbulent-boundary-layer 
interaction.  The  vortices  that  are  expected  to  appear  in  this  type  of  flow  are  visualized  in  the 
3-D  space  by  the  contours  of  the  eigenvalues  of  the  velocity  gradient  held.  In  addition,  three 
cross-sections  have  been  drawn  on  which  the  density  contours  are  displayed  (visualization 
of  shock  waves).  It  is  observed  that,  as  expected,  the  how  is  dominated  by  a  large  vortical 
structure,  which  lies  on  the  hat  plate  and  whose  core  has  a  remarkable  conical  shape  with  a 
battened  elliptical  cross-section.  Also  on  the  hat  plate,  on  the  side  of  the  main  vortex,  a  thin 
vortex  has  developed  in  the  direction  of  the  how.  This  is  not  an  independent  vortex,  but  the 
core  of  the  vorticity  sheet  which  lifts-off  the  surface,  along  the  separation  line,  and  rolls  up 
to  form  the  conical  vortex.  Along  the  vertical  hn  and  close  to  the  corner,  the  longitudinal 
vortex,  mentioned  by  Zheltovodov  [152,  153]  and  Kubota  and  Stollery  [146],  is  seen.  It  also 
develops  quasi-conically,  but  with  a  smaller  rate  of  increase,  compared  to  the  prime  conical 
structure.  In  the  lower  part  of  the  hgure,  a  cross-section  of  the  vortices  is  shown.  There, 
it  is  indicated  that  indeed  the  hat  ground  vortex  constitutes  the  initial  part  of  the  primary 
conical  vortex. 

In  the  conical  projections  of  Alvi  and  Settles  [149]  the  how  seems  to  be  quasi-conical 
rather  than  conical.  Then  the  question  arises,  what  causes  this  deviation?  In  this  context 
Panaras  [150]  observes  that  the  different  rate  of  thickening  of  the  conical  vortex  and  of  the 
boundary  layer  of  the  plate  is  expected  to  affect  the  conical  similarity  adversely.  For  studying 
this  effect  Panaras  [150]  projected  conically  the  sections  (i)  and  (ii)  shown  in  Fig.  58  on  the 
outflow  section  (iii).  The  iso-Mach  lines  have  been  used  as  visualization  parameters.  The 
results  of  that  correlation  are  shown  in  Fig.  59.  It  is  observed  in  Fig.  59b,  where  section  (i)  is 
conically  projected  on  the  outhow  section,  that  while  good  coincidence  is  observed  between 
the  separation  bubbles,  the  swept  shocks  and  the  shock  triple-point,  the  feet  of  the  A-shocks 
are  fairly  correlated,  especially  the  separation  shocks,  while  there  is  no  equivalence  at  all 
between  the  boundary  layers  in  the  two  cross-sections.  More  particularly,  the  boundary 
layer  of  the  hrst  cross-section  is  about  40%  thicker  than  the  boundary  layer  of  the  third 
section.  This  large  difference  in  the  scaling  is  due  to  the  lower  rate  of  development  of  the 
boundary  layer,  compared  to  that  of  the  vortex. 

The  deviation  from  the  conical  behaviour  is  smaller  if  the  two  cross-sections  that  are 
conically  correlated  are  closer.  This  is  demonstrated  in  Fig.  59c,  where  section  (ii)  is  conically 
projected  on  section  (iii).  Panaras  [150],  has  also  presented  quantitative  data  which  verify 
the  observation  done  in  Fig.  59b  that  the  deviation  from  conical  similarity  is  greater  at  the 
part  of  the  how  between  the  separation  shock  and  the  plate. 

In  addition  to  the  development  of  the  how  model  depicted  in  Fig.  57,  Alvi  and  Settles 
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[148]  have  also  discovered  some  peculiar  features  of  the  phenomenon  of  secondary  sepa¬ 
ration  in  a  single  fin  flow.  Indeed,  while  in  a  strong  swept  shock  laminar  boundary- layer 
interaction,  a  secondary  reattachment  line  appears  between  the  primary  separation  and  reat¬ 
tachment  lines,  their  experiments  indicated  that  no  secondary  reattachment  line  is  visible  in 
a  turbulent  interaction.  In  addition  to  this  difference,  between  laminar  and  turbulent  flows, 
Zheltovodov  [152,  153,  154],  Zubin  and  Ostapenko  [147]  and  Alvi  and  Settles  [148]  report 
another  one,  which  perhaps  is  more  fundamental.  While  in  laminar  flows  the  secondary 
separation  phenomenon  appears  progressively  with  increasing  shock  strength,  in  the  turbu¬ 
lent  fin/plate  flows  the  existence  of  a  secondary  separation  line  in  the  surface  flow  pattern 
has  been  detected  only  in  interactions  of  moderate  strength  and  it  disappears  in  stronger 
ones.  A  typical  skin-friction  pattern  for  a  single  fin  flow  with  secondary  separation  is  shown 
in  Fig.  60,  after  Panaras  and  Stanewsky  [155].  In  this  figure  it  is  demonstrated  that  the 
secondary  vortex  is  gradually  formed  along  the  secondary  separation  line.  No  secondary 
reattachment  line  is  observed.  We  note  in  Fig.  60  that  in  the  surface  flow  pattern  the  sepa¬ 
ration  and  the  reattachment  lines  are  curved  and  not  straight  close  to  the  apex.  The  initial 
region  of  deviation  from  conical  behavior  has  been  called  the  ”  inception  zone” . 

A  reasonable  explanation  of  the  strange  behavior  of  the  secondary  separation  in  a  single 
fin  flow  has  been  given  by  Zheltovodov  et  al  [154,  156].  According  to  them,  the  secondary 
separation  line  first  appears  once  the  interaction  has  achieved  a  certain  strength,  showing  up 
in  the  conical  region  but  not  in  the  inception  zone.  Its  spanwise  extent  grows  with  increasing 
shock  strength  but  then  diminishes  again,  eventually  appearing  only  in  the  inception  zone 
and  then  disappearing  altogether.  Secondary  separation  then  reappears  in  the  strongest 
interaction  observed  to  date  (Moo  =  4,  a  =  30.6°),  but  in  a  different  position,  closer  to  the 
fin  than  previously.  The  experimental  results  of  Zheltovodov  et  al  [156]  demonstrate  that 
the  initial  behaviour  of  secondary  separation  is  related  to  laminar,  transitional  and  then 
turbulent  reverse  flow  in  the  swept  separation  bubble.  They  also  ascribe  the  reappearance  of 
secondary  separation  to  the  development  of  supersonic  turbulent  reverse  flow  in  the  separated 
region  with  an  imbedded  normal  shock  wave  of  critical  strength  (with  the  pressure  ratio 
P2/P1  >  1.5  —  1.6)  for  turbulent  flows  to  specify  the  conclusion  by  Zubin  and  Ostapenko 
[147]  regarding  the  reason  of  disappearance  and  the  condition  of  reappearance  of  secondary 
separation.  Regarding  the  secondary  reattachment  line,  Zheltovodov  et  al  [156]  obtained 
evidence  of  its  existence  only  in  the  aforementioned  extremely  strong  interaction  (very  close 
to  the  secondary  separation  line).  The  work  of  Zheltovodov’s  group  is  summarized  in  Fig.  61. 
They  have  identified  six  regimes  for  the  3-D  single  £11  flowfield  depending  on  the  strength  of 
the  shock  wave.  In  Regime  I,  the  boundary  layer  is  unseparated  and  no  convergence  of  surface 
streamlines  is  observed.  In  Regime  II,  the  streamlines  turn  approximately  parallel  to  the 
inviscid  shock  but  do  not  form  a  line  of  coalescence.  I11  Regime  III,  a  primary  separation  line 
(SI)  forms  corresponding  to  the  coalescence  of  the  surface  streamlines.  An  attachment  line 
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(Rl)  forms  near  the  fin-plate  junction.  A  line  of  secondary  separation  (S2)  appears  located 
between  the  primary  separation  and  attachment  lines.  In  Regimes  IV  to  VI,  the  secondary 
separation  (S2)  and  attachment  lines  (R2)  disappear  and  finally  reappear.  Using  sand-grain 
roughness  mounted  along  the  primary  attachment  line  Rl,  Zheltovodov  et  al  [156]  showed 
experimentally  that  the  secondary  separation  can  be  suppressed  by  such  way  in  the  regimes 
III  and  IV,  but  not  in  the  regime  VI.  These  experiments  support  the  conclusion  regarding 
the  critical  strength  of  the  imbedded  normal  shock  wave  which  causes  the  reappearance  of 
secondary  turbulent  separation  in  the  regime  VI. 

Though  the  numerical  simulations  have  considerably  contributed  to  the  verification  of 
the  early  hypotheses  regarding  the  structure  of  the  single  fin  flows,  the  accuracy  of  the 
predictions  was  not  exceptional  in  the  early  90’s.  The  capability  of  numerical  simulation  of  3- 
D  shock- wave /turbulent  boundary-layer  interactions  has  been  assessed  by  Knight  [157],  who, 
to  that  purpose,  examined  five  specific  configurations  (sharp  fin,  blunt  fin,  cylinder /flare, 
swept  compression  corner  and  crossing  shocks)  at  Mach  numbers  from  two  to  eight.  Knight 
[157]  concluded  that  the  pitot  pressure,  yaw  angle  and  surface  pressure  are  predictable  with 
reasonable  accuracy  using  algebraic  or  two-equation  turbulence  models,  however  the  surface 
heat  transfer  is  not  accurately  predicted  in  strong  interactions.  In  the  particular  case  of 
the  sharp- fin/plate  configuration,  Kim  et  al  [158]  have  performed  a  joint  experimental  and 
computational  study  of  skin  friction  in  weak-to-strong  interactions  at  =  3  and  4.  In 
their  Navier-Stokes  calculations  they  tested  algebraic  turbulence  models  and  the  k—e  model, 
integrated  to  the  wall  or  employing  the  wall-function  technique.  They  have  found  that  their 
computations  agree  well  with  the  data  for  moderate  interaction  strengths,  but  systematically 
underpredict  the  data  with  increasing  interaction  strength.  Also,  the  secondary  separation 
line,  which,  according  to  the  experimental  data,  exists  in  the  surface  flow-pattern  has  not 
been  predicted  in  their  calculations. 

In  a  recent  work  Panaras  [2]  studied  the  structure  of  the  separation  vortex  in  a  strong 
swept  shock- wave/turbulent  boundary- layer  interaction,  for  the  purpose  of  explaining  why 
current  turbulence  models  fail  to  predict  accurately  this  type  of  flow.  One  of  the  sharp- 
fin/plate  flows  examined  by  Kim  et  al  [158]  was  used  as  model.  After  the  validation  of  the 
results,  by  comparison  with  appropriate  experimental  data,  he  studied  the  flowheld  by  means 
of  stream  surfaces  which  start  at  the  inflow  plane,  within  the  undisturbed  boundary  layer, 
and  which  are  initially  parallel  to  the  plate.  Each  of  these  surfaces  has  been  represented  by 
a  number  of  streamlines. 

Calculation  of  the  spatial  evolution  of  some  selected  stream  surfaces  has  revealed  that 
the  inner  layers  of  the  undisturbed  boundary  layer,  where  the  eddy  viscosity  is  high,  wind 
around  the  core  of  the  vortex.  However,  the  outer  layers,  which  have  low  turbulence,  rotate 
over  the  vortex  and  penetrate  into  the  separation  bubble  at  the  reattachment  region  forming 
a  low  turbulence  tongue,  which  lies  along  the  plate,  underneath  the  vortex  (Fig.  62).  The 
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intermittency  of  the  air  that  constitutes  the  tongue  and  the  outer  layers  of  the  vortex,  is 
very  small,  i.e.,  the  flow  is  almost  laminar  there.  At  the  initial  stage  of  development  the 
conical  vortex  is  completely  composed  of  turbulent  air,  but  gradually  downstream,  as  it 
grows  linearly,  the  low  turbulence  tongue  is  formed.  Panaras  [2]  observed  that  an  increase 
of  the  strength  of  the  interaction  results  in  the  folding  of  higher  layers  around  the  vortex, 
which  are  purely  inviscid.  For  example,  in  the  very  strong  interaction  shown  in  Fig.  63a,  the 
external  layers  of  the  separation  vortex  and  a  major  part  of  the  tongue  are  purely  inviscid.  In 
the  other  extreme,  in  a  weak  interaction  no  low-turbulence  tongue  is  formed  (Fig.  63b).  The 
existence  of  the  low-turbulence  tongue  underneath  the  conical  separation  vortex  of  a  strong 
swept  shock/turbulcnt  boundary  layer  interaction  creates  a  mixed-type  separation  bubble: 
turbulent  in  the  vicinity  of  the  separation  line  and  almost  laminar  between  the  reattachment 
line  and  the  secondary  separation  vortex.  This  type  of  separation  is  difficult  to  be  simulated 
accurately,  since  current  turbulence  models  are  based  on  the  physics  of  2-D  flows,  where  in 
a  separation  bubble  the  whole  recirculation  region  is  turbulent. 

For  improving  the  accuracy  of  the  numerical  predictions  of  swept  shock- wave/turbulent 
boundary-layer  interactions,  Panaras  [159]  considered  the  aforementioned  physical  feature 
in  the  derivation  of  a  new  equation  for  the  calculation  of  the  eddy- viscosity  coefficient  in  the 
region  of  the  separation  vortex.  He  followed  the  Baldwin-Lomax  [151]  formulation,  because 
it  is  easily  implemented  in  a  Navier-Stokes  computation  scheme.  The  modified  turbulence 
model  was  used  for  the  computation  of  some  of  the  test  cases  compiled  by  Settles  and  Dod¬ 
son  [8].  The  agreement  with  the  experimental  evidence  was  very  good.  In  addition,  the 
results  were  given  to  Knight  and  Degrez  [1]  who  compared  them  with  similar  calculations, 
contributed  by  other  researchers,  employing  different  turbulence  models,  in  order  to  assess 
the  capability  for  numerical  simulation  of  3-D  shock  wave  turbulent  boundary  layer  interac¬ 
tions.  The  turbulence  modeling  proposed  by  Panaras  [159]  was  the  only  one  that  predicted 
secondary  separation. 

In  principle,  the  eddy  viscosity  equation  developed  by  Panaras  [159]  is  appropriate  for 
applications  to  other  3-D  flows  characterized  by  the  appearance  of  extensive  crossflow  sep¬ 
aration  (generation  of  crossflow  vortices),  like  in  subsonic  or  supersonic  flows  about  slender 
bodies  or  delta  wings  at  high  incidence.  Very  good  agreement  with  the  experiments  was 
observed  for  the  Mach  2.0  flow  about  an  ogive-cylinder  at  a  =  10°  and  20°,  and  the  subsonic 
flow  about  a  6:1  ellipsoid  at  incidence  a  =  30°  as  described  in  Panaras  [160]. 

Very  recently,  Tliivet  et  al  [161]  in  an  effort  to  reduce  the  high  level  of  turbulent  kinetic 
energy  observed  within  the  separation  vortex  of  a  swept  shock  turbulent  boundary- layer 
interaction,  included  in  a  standard  k  —  cu  model  the  non-linear  correction  of  Durbin  [162], 
According  to  this  correction  in  regions  of  high  strain  the  eddy  viscosity  is  reduced.  Tliivet  et 
al  [161]  found  that  this  correction  improves  considerably  the  results  of  simulation  of  single 
fin  flows.  Their  results  are  included  in  this  review  and  are  presented  in  the  next  section. 
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Herein  we  just  note  that  according  to  their  results,  the  turbulent  kinetic  energy  within 
the  separation  vortex  is  considerably  smaller  than  that  given  by  the  standard  k  —  u:  model, 
resulting  in  a  vortex  of  larger  size.  It  is  remarkable  that  near  the  wall,  between  the  primary 
attachment  and  the  secondary  separation,  the  turbulent  energy  is  almost  zero,  exactly  as 
it  should  be  according  to  the  flow  model  of  Panaras  [2],  Indeed,  comparison  of  the  flow 
model  shown  in  Fig.  63a  and  the  cross-section  of  the  turbulent  kinetic  energy  calculated  by 
Thivet  (Fig.76)  reveals  a  remarkable  similarity.  In  both  cases  almost  laminar  flow  exists 
near  the  wall  underneath  the  separation  vortex,  and  in  both  cases  the  separation  vortex, 
outside  of  the  core,  consists  of  air  of  small  intermittency  or  of  small  turbulence,  respectively. 
Their  only  difference  is  the  core  itself.  In  the  model  of  Panaras,  the  core  includes  the  inner 
turbulent  part  of  the  incoming  boundary  layer,  while  in  Thivet’s  results  the  core  has  almost 
zero  turbulent  kinetic  energy.  However,  since  the  core  occupies  a  small  part  of  the  separation 
vortex,  its  effect  on  the  overall  flowfield  is  small.  In  conclusion,  the  comparison  of  Figs.  63a 
and  76  supports  the  view  that  actually  the  weakly  nonlinear  Wilcox-Durbin  model,  applied 
by  Thivet  et  al  [161],  predicts  a  flowfield  that  basically  is  similar  to  that  proposed  by  Panaras 
[2]- 

4.1  Cases  1-3:  Mach  3,  a  =  15°  and  Mach  4,  a  =  20°  and  30.6° 

Thivet  [161]  simulated  three  different  configurations  of  the  3-D  single  fin.  The  flow  conditions 
are  listed  in  Table  14.  Experimental  data  of  Zheltovodov  et  al  [143,  154]  and  Kim  et  al  [158] 
is  available  for  Cases  1,3  and  2,  respectively.  The  computations  were  performed  at  the  same 
flow  conditions  as  the  experiment.  The  inviscid  pressure  ratio  p^jpi  is  2.82,  5.21  and  9.5  for 
Cases  1  to  3,  respectively. 

Details  of  the  computation  are  presented  in  Table  15.  Two  different  turbulence  models 
were  examined,  namely,  the  linear  Wilcox  k  —  uj  model  [163]  (WI)  and  a  weakly  nonlinear 
Wilcox  k  —  to  model  [164]  (WD+).  The  WD+  model  is  an  extension  to  compressible  flows 
of  the  nonlinear  correction  of  Durbin  [162],  I11  the  WD+  model,  the  coefficient  in  the 
turbulent  eddy  viscosity  pt  =  pc^k/uj  is  defined  as  =  min(c°,  y/c°/s)  where  s  =  S/u, 

S  =  sj2SijSij  —  § SkkSkk  and  Stj  =  §( dUi/dxj  +  dUj/dxi).  Computations  were  performed 
using  the  GASPex  code  [165].  Inviscid  fluxes  are  computed  to  third  order  using  the  Roe 
scheme  with  MUSCL  reconstruction  and  a  min-mod  limiter.  Viscous  and  heat  transfer  terms 
are  computed  using  second  order  central  differences.  A  typical  grid  is  shown  in  Fig.  64. 

The  surface  streamlines  are  displayed  in  Figs.  65  (WI),  66  (experiment)  and  67  (WD+) 
for  Case  1,  and  in  Figs.  68  (WI),  69  (experiment)  and  70  (WD+)  for  Case  3.  The  angle  of 
the  primary  separation  line  Si  is  underestimated  by  the  linear  (WI)  model  by  1.4°  to  4°, 
and  overestimated  by  the  weakly  nonlinear  (WD+)  model  by  1.4°  to  2.5°.  The  angle  of  the 
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primary  attachment  line  Ri  is  accurately  predicted  by  both  models.  In  Case  1  (Regime  III), 
the  weak  secondary  separation  S2  observed  in  the  experiment  (Fig.  66)  is  absent  in  the  com¬ 
putations  with  WI  model  (Fig.  65).  For  the  WD+  model  (Fig.  67),  the  distinct  changes 
in  the  curvature  of  surface  streamlines  occurring  midway  between  the  hn  and  primary  sep¬ 
aration  line  are  similar  to  incipient  secondary  separation  conditions  which  occurred  in  the 
experiment  at  a  lower  £11  angle.  I11  Case  3  (Regime  VI),  the  strong  secondary  separation  S2 
and  attachment  R2  lines  observed  in  the  experiment  (Fig.  69)  are  evident  in  the  WD+  model 
(Fig.  70)  but  are  entirely  absent  in  the  WI  model  (Fig.  68).  This  is  a  significant  achievement 
of  the  WD+  model. 

The  computed  and  experimental  pressure  distributions  are  displayed  in  Figs.  71  to  73. 
The  WD+  model  shows  a  significant  improvement  compared  to  the  WI  model,  and,  in 
particular,  accurately  predicts  the  appearance  of  the  secondary  shock  wave  associated  with 
the  secondary  separation.  As  is  seen  from  Figs.  70  and  and  73,  the  reappearance  of  the 
secondary  separation  in  the  regime  VI  has  been  achieved  in  the  computations  with  WD+ 
model  when  the  strength  of  the  imbedded  normal  shock  wave  was  P2/P1  =  2,  as  in  the 
experiment.  This  value  is  higher  compared  to  the  critical  value  P2/P1  =  1-5  —  1.6  which  is 
necessary  to  cause  the  appearance  of  turbulent  separation. 

The  computed  and  experimental  skin  friction  coefficient  for  Case  2  is  shown  in  Fig.  74. 
The  peak  skin  friction  for  the  WD+  model  is  reduced  by  35%  compared  to  the  WI  model. 
Two  series  of  measurements  are  shown,  where  the  solid  symbols  •  correspond  to  corrected 
values  of  the  skin  friction  [166].  The  WD+  model  displays  significantly  better  agreement 
with  experiment  than  the  WI  model.  In  accordance  with  the  surface  skin  friction  prediction 
(Fig.  74)  as  well  as  the  surface  pressure  distribution  considered  by  Tliivet  et  al  [161]  for 
regime  IV  (M  =  4,  a  =  30.6°),  secondary  separation  didn’t  exist  in  the  computations,  as  in 
experiments,  when  the  strength  of  the  embedded  shock  wave  was  smaller  that  the  critical 
value  for  the  turbulent  separation 

The  improvement  achieved  by  the  WD+  model  is  associated  with  a  significant  reduction 
in  the  peak  k  within  the  primary  separation  vortex  core  and  in  the  near  wall  flow  penetrating 
to  the  place  of  formation  of  the  secondary  separation  line.  Figs.  75  and  76  display  computed 
turbulence  kinetic  energy  contours  for  the  WI  and  WD+  models,  respectively.  The  WD+ 
model  reduces  the  peak  k  within  the  core  by  more  than  a  factor  of  two,  thereby  increasing 
the  size  of  the  interaction.  These  results  also  support  the  conclusion  of  Zhcltovodov  et  al 
[143,  156,  154]  that  the  turbulence  amplification  in  the  reversed  flow  may  be  the  reason  for 
the  disappearance  of  the  secondary  separation  in  the  regime  IV  -  V  (Fig.  61). 

I11  summary,  the  following  conclusions  can  be  made: 

•  The  weakly  nonlinear  Wilcox-Durbin  (WD+)  model  provides  a  significant  improvement 
in  prediction  of  the  3-D  single  fin  interaction  in  comparison  with  the  Wilcox  (WI) 
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model.  In  particular, 

—  The  secondary  separation/attachment  and  secondary  shock  are  captured  for  the 
Mach  4,  a  =  30.6°  configuration  (Regime  VI).  The  linear  Wilcox  model  does  not 
predict  secondary  separation/attachment  and  the  secondary  shock  for  this  case. 

—  The  surface  pressure  distribution  is  in  close  agreement  with  experiment  for  all 
three  cases,  except  for  the  position  of  upstream  influence  (he.,  the  location  of 
initial  pressure  rise)  for  the  strongest  interaction. 

•  The  improvement  in  prediction  by  the  Wilcox-Durbin  model  is  associated  with  a  re¬ 
duction  in  the  computed  turbulence  kinetic  energy  k  in  the  separation  vortex. 

5  3-D  Double  Fin 

The  3-D  double  hn  interaction  (also  known  as  the  crossing  shock  wave  /  turbulent  boundary 
layer  interaction  )  shown  in  Fig.  5  is  a  canonical  validation  case  for  CFD  modeling  of  3- 
D  shock  wave  turbulent  boundary  layer  interactions.  An  equilibrium  supersonic  turbulent 
boundary  layer  approaches  a  channel  defined  by  two  semi-infinite  fins  of  angles  aq  and  a2 
separated  by  a  distance  We  at  the  entrance.  Additional  geometric  parameters  include  the 
minimum  (throat)  width  Wt  and  the  offset  (if  any)  of  the  midpoint  of  the  throat  relative  to 
the  midpoint  of  the  entrance. 

The  first  numerical-experimental  studies  of  such  configuration  with  small  and  moderate 
fin  inclination  angles  were  conducted  more  than  a  decade  ago  by  Mee,  Stalker  and  Stollery 
[167].  Since  that  time,  interest  in  this  problem  has  greatly  increased  due  to  its  relationship  to 
scramjet  propulsion  since  the  crossing  shock  interaction  constitutes  a  geometrical  simplifica¬ 
tion  of  a  hypersonic  sidewall  compression  inlet.  Symmetric  interactions  (i.e.,  aq  =  a2)  have 
been  studied  by  [6,  159,  164,  168,  169,  170,  171,  172,  173,  174,  175,  176,  177,  178,  179,  180, 

181,  182,  183,  184,  185,  186,  187,  188,  189,  190,  191,  192,  193,  194,  195,  196,  197,  198,  199, 

200,  201,  202,  203,  204,  205,  206,  207,  208,  209].  Asymmetrical  interactions  (i.e.,  aq  ^  a2) 

have  been  studied  by  [6,  164,  167,  169,  171,  172,  173,  174,  179,  187,  192,  193,  196,  197,  198, 

199,  202,  203,  204,  208,  209,  210,  211,  212,  213,  214,  215,  216,  218],  These  studies  cover  the 
Mach  number  range  from  1.85  to  8.3  and  fin  angles  from  —3°  to  15°.  Special  emphasis  in 
most  papers  is  on  the  analysis  of  the  pressure  distribution  and  the  structure  of  the  limiting 
streamlines  on  the  plate  surface  (surface  flow  pattern)  in  the  interaction  region. 

The  most  complete  cycle  of  experimental  research  of  the  symmetrical  interaction  be¬ 
tween  crossing  shock  waves  with  a  turbulent  boundary  layer  for  M, ^  =  3  and  4  is  presented 
in  [172,  174,  177,  178,  184,  185,  186].  In  addition  to  the  patterns  of  the  limiting  streamlines 
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and  pressure  distribution  on  the  plate  surface,  these  papers  present  the  data  on  the  distri¬ 
bution  of  the  skin  friction  coefficient,  total  pressure  fields  and  local  angles  of  flow  direction, 
as  well  as  the  optical  visualization  of  the  3-D  flowfield.  Similarly,  fairly  complete  mea¬ 
surements  of  surface  pressure  distribution  and  heat  transfer  together  with  fields  of  various 
parameters  under  the  conditions  of  symmetrical  interaction  are  analyzed  for  two  situations 
(Mqo  —  8.3,  a  =  10°  and  15°)  in  [179,  180, 181, 182],  Additional  systematic  information  about 
the  surface  heat  transfer  together  with  surface  pressure  distribution  and  limiting  streamlines 
pattern  was  obtained  at  =  4  for  different  cases  of  symmetrical  and  asymmetrical  in¬ 
teractions  in  [187,  194,  196,  216].  The  pioneering  study  [219]  of  a  symmetrical  interaction 
for  Mqo  =  3.85  near  the  a  =  15°  fins  mounted  on  the  surface  of  a  2-D  compression  corner 
(a  =  10°)  is  also  noted.  The  data  on  the  gas  dynamic  structure  of  this  more  complex  flow 
obtained  on  the  basis  of  optical  visualization  and  the  limiting  streamlines  on  the  surface  are 
important  for  testing  the  up-to-date  calculations  and  searching  for  effective  configurations 
of  inlets  with  3-D  compression. 

Studies  dealing  with  unsteady  effects  are  infrequent,  but  nonetheless  extremely  impor¬ 
tant.  Batcho  et  al  [168]  were  apparently  the  first  to  measure  the  pressure  and  temperature 
fluctuations  on  the  surface  for  the  case  of  symmetrical  interaction.  These  measurements 
allowed  the  authors  to  assume  the  existence  of  a  relationship  between  the  fluctuations  of 
these  parameters  and  heat  transfer  intensity  fluctuations.  More  detailed  measurements  of 
pressure  fluctuations  were  performed  later  in  [169,  172], 

It  should  be  noted  that  fins  with  flat  internal  faces  that  form  a  converging  half-channel 
were  used  in  most  experiments,  except  for  [179,  180,  181,  182,  187,  194,  196,  216].  There 
was  only  one  intersection  of  shock  waves  because  of  the  threat  of  half-channel  blockage,  and 
the  downstream  region  of  investigation  on  the  plate  was  limited.  In  the  above  papers,  the 
side  faces  had  inflected  generatrices  with  a  constant-width  channel  behind  the  inflection  (see 
Fig.  5).  Specific  features  of  the  flow  in  these  conditions  are  determined  by  the  boundary 
layer  interaction  with  a  sequence  of  crossing  shocks  and  expansion  waves. 

Advanced  numerical  studies  of  these  flows  are  based  on  the  numerical  solution  of  the 
full  Reynolds  averaged  Navier-Stokes  equations  (RANS)  with  the  use  of  various  turbulence 
models:  algebraic  (Baldwin- Lomax  and  its  modifications)  and  differential  models  (k— e,  k— u>, 
k  —  l  and  their  modifications),  and  also  the  full  Reynolds  stress  equations.  The  calculations 
conducted  for  symmetrical  [159,  164,  170,  174,  175,  176,  177,  180,  181,  182,  183,  184,  186,  188, 
189,  190,  191,  192,  193,  199,  200,  201,  202,  203,  204,  205,  206,  208,  209,  210,  213,  214]  and 
asymmetrical  [164,  187,  192,  193,  197,  199,  202,  203,  204,  208,  209,  211,  212,  213,  214,  215, 
216,  218]  interactions  and  for  symmetrical  fins  mounted  on  the  compression  surface  [190,  219]. 
As  it  follows  from  review  of  investigations  of  the  flow  for  symmetrical  and  asymmetrical 
interactions,  the  advanced  RANS  calculations  with  various  turbulence  models  (algebraic 
or  two-parameter  differential  models,  and  also  Reynolds  stress  equations)  allow  a  correct 
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prediction  of  only  some  properties  of  these  flows.  For  example,  the  numerical  results  are 
in  good  agreement  with  experimental  data  for  the  shock-wave  structure  formed  in  the  half¬ 
channel  under  the  conditions  of  viscous-inviscid  interaction,  the  fields  of  various  parameters 
of  the  mean  flow,  the  pressure  distribution  on  the  surface  only  for  limited  range  of  the 
crossing  shock  wave  strength  (under  the  conditions  of  weak-  and  moderate-strength).  This 
agreement,  however,  significantly  decreases  as  the  strength  of  crossing  shock  waves  increases. 
The  accurate  prediction  of  surface  heat  transfer  and  skin  friction  remains  a  challenging 
problem. 

To  characterize  some  specific  features  of  the  crossing  shock  wave/turbulent  boundary 
layer  interaction  Fig.  77  shows  a  3-D  perspective  view  of  the  flowheld  model  with  schematic 
images  of  the  interaction  structure  at  some  specific  cross  sections  (I  to  III)  in  the  vicinity 
of  symmetrical  15°  x  15°  double  hn  configuration  at  M, ^  =  4  developed  by  Garrison  and 
Settles  [173,  177,  178]  on  a  basis  of  experimental  planar  laser  scattering  images.  Based  on  the 
inherent  symmetry  of  the  crossing-shock  interaction,  only  half  of  interaction  is  shown  in  the 
figure.  Cross  section  I  occurs  upstream  of  the  intersection  of  the  two  single  hn  interactions 
meet  and  thus  represents  two  separate  single  hn  interactions.  Incident  separation  (1),  rear 
(2)  and  “inviscid”  (3)  shocks  with  a  slip  line  (4)  penetrating  from  the  triple  point  to  a  hn 
side/plate  cross  line  as  well  as  separation  vortex  (5)  under  the  bifurcated  shock  system  are 
visible  in  this  cross  section.  Cross  section  II  occurs  downstream  of  the  intersection  of  the  two 
separation  shocks.  To  understand  the  intersection  of  the  two  separate  single  hn  interactions, 
the  vertical  plane  of  symmetry  is  considered  an  inviscid  reflection  plane.  For  the  symmetrical 
crossing  shock  wave  interaction,  shock  waves  that  intersect  this  plane  must  rehect  from  it  to 
satisfy  continuity.  It  is  evident  that  in  cross  section  II  the  incident  separation  shock  reflection 
from  the  symmetry  plane  is  an  irregular  (he.,  Mach)  reflection.  The  Mach  reflection  results 
in  a  straight  shock  segment,  the  Mach  stem  (7),  which  spans  the  interaction  centerline,  a 
rehected  portion  of  the  separation  shock  (6),  and  newly  formed  triple  point  (10).  As  it  is 
seen  from  the  perspective  view  and  cross  section  III,  the  entire  incident  A-shock  structure 
rehects  from  the  center  plane  in  an  irregular  manner  and  remains  intact,  though  somewhat 
distorted,  propagating  away  from  the  center  toward  the  hn  surface.  The  additional  cross 
sections  that  characterize  details  of  the  complex  shock  crossings  between  the  cross  sections 
I  and  III  are  considered  in  [173,  177,  178].  As  shown  in  cross  section  III,  two  additional 
shock  segments  and  two  triple  points  are  observed  to  form  a  result  of  this  crossing.  One 
centerline  segment  (12)  spans  the  interaction  centerline  between  the  rehected  inviscid  shock 
waves  (9),  and  the  second  “bridge”  segment  (13),  connects  the  rehected  inviscid  shock  (9)  to 
the  rehected  separation  (6)  and  rear  (8)  shocks.  A  “mushroom-shaped”  separation  structure 
(5)  is  formed  in  the  vicinity  of  the  centerline  in  the  last  cross  section  from  the  twin  singlc-hn 
separations  which  begin  the  crossing-shock  interaction. 

An  example  of  the  howheld  for  a  M, ^  =  4,  15°  x  15°  symmetric  double  hn  interaction  is 
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shown  in  Fig.  78  by  the  stream  ribbons  together  with  the  footprint  of  each  regime  (coherent 
feature)  (Fig.  79)  in  accordance  with  RANS  computations  by  Gaitonde  and  Shang  [189]  and 
Gaitonde  et  al  [192,  193,  197]  performed  for  the  test  configuration  investigated  experimen¬ 
tally  by  Garrison  and  Settles  [178]  and  Zheltovodov  et  al  [187,  194],  The  turbulence  model 
employed  to  derive  the  eddy  viscosity  is  based  on  the  k  —  e  equations  [220,  221]  with  low 
Reynolds  number  terms  and  incorporates  a  compressibility  correction.  As  pointed  out  by 
Gaitonde  and  Shang,  different  coherent  features,  or  regimes,  can  be  identified  in  the  flow. 
For  example,  it  is  evident  from  Fig.  78  [189,  193,  197,  218]  that  the  incoming  boundary  layer 
separates  at  the  line  of  coalescence  SI  and  its  symmetric  counterpart  S 2  (not  marked).  This 
regime  does  not  reattach  thus  setting  up  an  open  structure  with  no  “trapped”  fluid.  Flow 
attaching  near  the  fin-plate  corners  at  Rl  and  its  symmetric  counterpart  R2  (not  marked) 
(see  also  Fig.  79)  can  be  categorized  into  several  regimes,  the  first  of  which  is  a  Vortex  In¬ 
teraction  (VI),  characterized  by  off-surface  stagnation  points  [192,  200,  201,  206]  and  which 
separates  from  the  downstream  side  of  51.  Fluid  attaching  further  downstream  forms  the 
entrainment  flow  (EF)  which  brings  high  speed  fluid  near  the  plate.  This  regime  separates 
prior  to  reaching  the  symmetry  plane,  and  two  Centerline  Vortices  (CV)  are  formed. 

The  several  lines  of  coalescence  and  divergence  observed  in  Fig.  79  have  direct  counter¬ 
parts  in  experimental  observation,  as  documented  previously  [173,  177,  178,  187,  194,  196, 
205,  206].  For  example,  in  accordance  with  experiments  by  Zheltovodov  et  al  [187,  194,  196] 
the  primary  separation  lines  Si  and  S2,  the  secondary  separation  lines  S3  and  54  and  the 
primary  attachment  lines  R\  and  R2  are  specific  for  15°  x  15°  double-fin  case  at  M ^  =  4 
(Fig.  80a).  The  distinct  separation  region  is  formed  in  the  throat  behind  the  centerline  sin¬ 
gular  saddle  point  Si  and  bounded  by  the  separation  lines  S$  and  Sq.  The  centerline  node 
point  Nl  was  observed  further  downstream.  The  line  S3  is  very  close  to  merge  with  S$  (and 
so  is  54  with  S6  on  the  opposite  side)  approximately  in  the  middle  of  the  central  separated 
zone. 

The  topology  of  surface  flow  pattern  computed  by  Gaitonde  and  Shang  on  a  basis  of 
the  k  —  e  turbulence  model  [192,  205,  206]  is  somewhat  different  (Fig.  80b, c).  The  primary 
separation  lines  intersect  at  the  centerline  node  N%,  and  therefore  no  fluidic  throat  forms  but 
in  the  experiment  another  centerline  node  Nl  is  observed.  (This  is  more  clearly  shown  in  an 
enlarged  fragment  in  Fig.  80c).  The  unsteady  intermittent  nature  of  the  experimental  flow 
upstream  of  the  singular  separation  point  Si  [194,  217],  not  predicted  by  computations,  may 
explain  this  discrepancy.  Additional  features,  which  were  not  observed  in  the  experiment, 
are  two  saddle  points  S±  and  S2  located  symmetrically  about  the  centerline  (Fig.  80c).  As 
shown  by  Schulcin  and  Zheltovodov  [195]  these  symmetric  saddle  points  appear  as  a  result 
of  the  complex  evolution  of  the  flow  at  higher  shock  waves  strength. 

The  calculations  do  not  reveal  the  signs  of  the  secondary  separation  lines  S3  and  54,  and 
the  central  separated  zone  is  more  compressed  by  the  flows  propagating  from  the  fins  and 
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noticeable  smaller  in  width  comparing  with  experiment  (Fig.  80a, b).  As  result,  the  computed 
surface  pressure  is  noticeably  overestimated  in  the  vicinity  of  the  channel  centerline  for 
X  =  x/d  >  17  (Fig.  80d).  As  was  concluded  in  [198,  205,  206],  this  discrepancy  may  be 
related  to  a  difference  in  the  turbulence  level  in  the  near  wall  region  of  the  secondary  flows 
between  the  primary  attachment  ( Ri  and  /?,2)  and  separation  (Si  and  S2)  lines  in  experiment 
and  computation. 

This  conclusion  is  supported  by  the  computations  of  Panaras  [159]  (Fig.  81b)  who  has 
used  his  modification  of  the  Baldwin-Lomax  turbulence  model  to  predict  a  similar  M <*>  =  4, 
15°  x  15°  crossing  shock  wave  turbulent  boundary  layer  interaction  investigated  by  Garrison 
and  Settles  [177,  178,  184,  186]  (Fig.  81a).  The  resultant  pattern  of  the  limiting  streamlines 
(Fig.  81b)  is  a  good  reproduction  of  the  specific  features  of  the  behavior  of  the  secondary 
separation  lines  S3  and  S4  under  the  conditions  of  the  laminar  secondary  flow,  the  scale  of 
the  central  separation  region,  the  centerline  pressure  distribution  (Fig.  81c)  and  skin  friction 
distribution  (Fig.  81c),  which  were  observed  by  Garrison  and  Settles  in  their  experiments. 
At  the  same  time,  the  computed  surface  flow  pattern  is  associated  with  the  node  in  the  apex 
of  the  centerline  separation  zone  and  two  symmetrical  saddle  points  as  in  the  computations 
by  Gaitonde  and  Shang  considered  above  which  was  not  observed  in  the  experiment  at  such 
conditions. 

Gaitonde  et  al  [192,  193,  197,  218]  performed  RANS  computations  with  the  k  —  e  model 
and  described  the  evolution  of  the  flow  model  with  interaction  strength  and  asymmetry.  An 
attempt  was  made  in  [193,  201,  218]  to  inhibit  the  transition  of  the  sidewall  boundary  layer 
by  limiting  the  production  of  turbulent  kinetic  energy  in  a  k  —  e  model.  It  resulted  in  a 
rectification  of  the  discrepancies  between  the  flowfield  visualization  and  simulations  in  the 
sidewall/shock- vortex  interaction  region  but  no  significant  effect  on  the  bottom  wall. 

Recent  studies  [159,  192,  193,  197,  199,  206,  211,  216]  distinctly  demonstrate  obvious 
limitations  of  various  turbulence  models  to  correctly  predict  different  properties  of  symmetric 
and  asymmetric  crossing  shock  wave  turbulent  boundary  layer  interaction  for  a  wide  range 
of  crossing  shock  strength.  They  have  stimulated  new  attempts  to  analyze  the  possibilities 
of  two-equation  k—tu  and  k—e  turbulence  models  and  some  of  their  modifications,  as  well  as 
modifications  of  the  algebraic  Baldwin-Lomax  model  by  Panaras  for  predicting  the  properties 
of  such  flows  over  a  wide  range  of  crossing  shock  wave  strength  [6,  164,  200,  201,  202,  203, 
204,  205,  208,  209,  210,  218], 
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5.1  Case  1:  Mach  4,  aq  =  15°,  aq  =  15° 

Thivet[164,  202,  204,  209]  simulated  two  configurations  of  the  3-D  double  fin6.  Experimental 
data  of  Zhcltovodov  et  al  [187,  194,  196,  216]  is  available  for  both  cases.  The  flow  conditions 
are  listed  in  Table  16.  The  computations  were  performed  at  the  same  flow  conditions  as 
the  experiment.  The  experimental  configuration  for  Case  1  is  shown  in  Fig.  5.  The  fins 
were  100  mm  high  and  192  mm  long.  Their  leading  edges  were  located  at  a  distance  of 
210  mm  from  a  leading  edge  of  the  plate.  The  boundary  layer  thickness  upstream  of  the 
channel  entrance  d  =  3.5  mm.  The  entrance  width  We  =  79.1  mm  and  the  throat  width 
Wt  =  32  mm. 

Details  of  the  computations7  are  presented  in  Table  17.  Five  different  turbulence  models 
were  considered,  namely,  the  (linear)  Wilcox  k  —  to  (WI)  model  [163],  and  four  variants  of 
the  Wilcox  model  based  on  the  requirement  for  weak  realizability [222]  of  the  computed 
Reynolds  stresses8.  The  variants  are  the  Wilcox-Moore  (WM),  Wilcox-Durbin  (WD)  and 
two  modifications  thereof  (WM+  and  WD+).  Details  of  the  models  are  presented  in  Tliivet 
et  al  [164],  A  detailed  grid  refinement  study  [164]  was  performed  to  quantify  the  uncertainty 
in  the  computations. 

The  surface  flow  pattern  topology,  skin  friction  lines  and  heat-transfer  coefficient  lines 
computed  by  Thivet  [164,  204,  209]  using  the  WI  and  WM+  models  are  presented  in  Fig.  82. 
The  computed  primary  separation  lines  ST  and  S 2  intersect  at  the  node  point  and  the 
fluidic  throat  does  not  exist  between  these  lines  (compared  with  the  experiment  where  the 
saddle  point  has  been  revealed  in  the  throat  between  these  lines,  see  Fig. 77).  The  secondary 
separation  lines  S3  and  S4  are  not  predicted  and  additional  separation  lines  S5  and  S 6, 
which  bound  the  central  separation  region,  appear  much  closer  to  the  centerline  than  in  the 
experiment.  Nevertheless,  the  width  of  this  separation  zone  is  wider  comparing  with  one  in 
the  computations  by  Gaitonde  and  Shang  (Fig.  80b, c). 

The  computed  and  experimental  surface  pressure  and  heat  transfer  coefficient  on  the 
Throat  Middle  Line9  (TML)  for  Case  1  are  shown  in  Fig.  83  for  the  WI  and  WM+  models. 
In  the  case  under  consideration,  two  computations  were  carried  out  with  WM+  model  on 
fine  and  very  fine  grids.  All  three  computations  with  different  models  underestimate  the 
location  of  the  pressure  rise  immediately  following  the  plateau  and  overestimate  the  peak 
pressure.  The  subsequent  expansion  (associated  with  the  corner  of  the  fins)  is  more  accu- 

6A  third  configuration  =  4,  cu  =  7°,a2  =  7°)  is  not  presented  here,  as  it  is  a  weaker  interaction.  It 
is  accurately  predicted  by  general  RANS  models. 

7Data  are  shown  for  the  finest  grid. 

8Realizability  implies  that  the  normal  kinematic  Reynolds  stresses  u'u',v'v'  and  w'w'  are  non-negative, 
and  the  Reynolds  shear  stresses  u'v'  ,u'w'  and  v'w'  are  bounded  by  the  Schwartz  inequality,  e.g.,  (u'v')2  < 
u'u'  v'v' . 

9The  Throat  Middle  Line  is  the  streamwise  datum  through  the  center  of  the  throat. 
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rately  predicted  by  the  WI  model.  The  realizability  correction  does  not  modify  the  computed 
topology  so  that  the  features  of  the  heat  transfer  distribution  remain  the  same.  Nevertheless, 
the  levels  are  significantly  modified,  essentially  because  of  the  lower  turbulence  levels  in  the 
WM+  solution,  which  allow  an  earlier  appearance  (he.,  farther  from  TMP)  of  the  separation 
lines  S§  and  Sq  than  in  the  WI  solution,  yielding  a  wider  central  separation  region  and  a 
lower  heat  transfer  between  the  maxima.  The  computed  peak  heat  transfer  significantly 
overestimates  the  experiment  for  all  models. 

In  summary,  the  following  conclusions  can  be  made: 

•  The  surface  pressure  is  predicted  with  reasonable  accuracy,  although  some  discrepan¬ 
cies  between  computed  and  experimental  surface  pressure  are  evident. 

•  The  peak  heat  transfer  is  not  accurately  predicted  by  any  of  the  models. 

5.2  Case  2:  Mach  4,  a\  =  7° ,0.2  =  11° 

An  example  of  the  mean  streamline  structure  is  shown  in  Fig.  84  from  Thivet  et  al  [164] 
for  a  Mach  4,  7°  x  11°  double  fin  interaction.  The  principal  features  include  the  separated 
boundary  layer,  vortex,  and  entrained  flow.  The  computed  surface  heat  transfer  coefficient 
contours  and  surface  streamlines  are  shown  in  Fig.  85  and  the  corresponding  experimental 
surface  streamline  pattern  in  Fig.  86.  The  primary  separation  lines  (denoted  S 1  and  S2  in 
Fig.  86)  are  evident  in  both  computation  and  experiment,  and  correspond  to  the  liftoff  of 
the  streamlines  from  the  surface  in  Fig.  84.  The  attachment  lines  (denoted  Ri  and  R2  in 
Fig.  86)  correspond  to  the  impingement  of  the  streamlines  which  demarcate  the  boundaries 
of  the  counter-rotating  vortex  pair  formed  by  the  separated  flow.  Secondary  separation  lines 
(S3  and  S4  in  Fig.  86)  and  attachment  line  (denoted  R3  in  Fig.  86)  are  also  evident. 

The  computed  and  experimental  surface  pressure  and  heat  transfer  coefficient10  on  the 
TML  for  Case  2  are  shown  in  Fig.  87.  The  surface  pressure  is  accurately  predicted  by  all 
models  with  only  slight  differences  in  the  plateau  pressure  at  x  =  80  — 100  mm.  Similar  good 
agreement  with  experiment  and  insensitivity  of  the  predictions  to  the  turbulence  model  are 
observed  in  the  spanwise  pressure  distributions  at  three  streamwise  locations  (Fig.  88).  This 
is  consistent  with  previous  results  [1]  which  showed  a  marked  insensitivity  of  the  computed 
surface  pressure  to  the  turbulence  model  employed. 

The  computed  surface  heat  transfer  coefficient  on  the  TML  (Fig.  87)  and  at  three  stream- 
wise  locations  (Fig.  89)  significantly  overpredicts  the  peak  experimental  value  for  all  models. 
This  is  similar  to  previous  results  using  other  RANS  models  which  did  not  incorporate  a 
specific  realizability  effect  [1], 

10The  heat  transfer  coefficient  =  qw / PooUooCp(Tw  —  Taw)  where  Taw  is  the  adiabatic  wall  temperature. 
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Thivet  examined  the  effect  of  imposing  the  turbulence  length  scale  limitation  proposed 
by  Coakley  and  Huang  [223].  The  turbulent  eddy  viscosity  =  pc^^i  and  the  turbulence 
length  scale  i  is  defined  by 


^  =  min 


4>g,  k1/2/uj 


where 


dog 


K 


where  k,  =  0.41  is  von  Karman’s  constant  and  d  is  the  distance  normal  to  the  wall.  The 
dissipation  term  in  the  transport  equation  for  k  is  modeled  using  to  =  kx!2i.  Fig.  90  displays 
the  ratio  of  the  computed  turbulence  length  scale  i  (without  the  Coakley  and  Huang  correc¬ 
tion)  to  the  logarithmic  expression  t']0g  for  Case  2  at  x  =  46  mm  for  six  spanwise  locations. 
The  transverse  location  zr  is  zero  at  the  fin  and  unity  at  the  location  corresponding  to  the 
inviscid  shock  generated  by  the  11°  fin.  It  is  evident  that  the  ratio  does  not  significantly 
exceed  unity,  thereby  implying  that  the  Coakley  and  Huang  correction  would  have  little 
effect  on  the  computed  solution. 

Thivet  also  examined  the  effect  of  incorporating  a  compressibility  correction  [163]  in  the 


k  —  u;  model11.  The  threshold  value  of  the  turbulence  Mach  number  Mt  =  \j2k/a  (where 
a  is  the  local  speed  of  sound)  is  0.25  for  the  compressibility  correction.  Below  this  value, 
the  correction  is  inoperative.  The  maximum  value  of  Mt  using  the  WI  model  is  0.3  yielding 
an  increase  of  only  4%  in  the  dissipation  term,  thereby  implying  that  this  correction  (if 
implemented)  would  have  negligible  effect  on  the  computation  using  the  WI  model.  The 
values  of  Mt  for  the  WD,  WD+  and  WM+  models  are  everywhere  below  the  threshold, 
and  therefore  the  compressibility  correction  (if  implemented)  would  have  no  effect  on  these 
results. 

In  summary,  the  following  conclusions  can  be  made: 


•  The  computed  surface  pressure  is  accurately  predicted  by  the  five  different  turbulence 
models  for  Case  1  (aq  =  7°,  «2  =  11°). 

•  The  peak  heat  transfer  is  not  accurately  predicted  by  any  of  the  models. 

•  The  proposed  length  scale  correction  of  Coakley  and  Huang  [223],  and  the  compress¬ 
ibility  modifications  of  Zeman  [224]  or  Sarkar  [225],  would  not  significantly  change  the 
computed  flowhelds. 

11The  compressibility  correction  concept,  introduced  by  Zeman  [224]  and  Sarkar  [225]  to  account  for  the 
observed  decrease  in  the  spreading  rate  of  compressible  mixing  layers,  leads  to  an  undesirable  reduction  in 
skin  friction  for  turbulent  boundary  layers  [163,  223,  226]. 
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5.3  Case  3:  Mach  5,  aq  =  18°,  aq  =  18° 

Schmisseur  and  Gaitonde  [200]  and  Panaras  [227]  simulated  a  Mach  5,  aq  =  «2  =  18°  double 
fin  at  Res  =  1.4  x  105  and  isothermal  wall  conditions  ( Tw/Taw  =  0.76).  Experimental  data  of 
Schiilcin  and  Zhcltovodov  [195,  207]  (an  additional  description  is  presented  in  Zheltovodov 
and  Maksimov  [228])  is  available  for  this  case.  The  flow  conditions  are  shown  in  Table  18. 
The  computations  were  performed  at  the  same  flow  conditions  as  the  experiment.  The 
experimental  configuration  is  shown  in  Fig.  91. 

Details  of  the  computations  are  presented  in  Table  19.  Schmisseur  and  Gaitonde  utilized 
the  k—e  model  [220,  221]  with  low  Reynolds  number  and  compressibility  corrections.  Panaras 
[2]  employed  a  modified  Baldwin-Lomax  model. 

The  experimental  surface  oil  flow  is  presented  in  Fig.  92,  and  the  computed  surface 
streamlines  in  Fig.  93  (Schmisseur  and  Gaitonde)  and  Fig.  94  (Panaras).  The  experiment 
displays  a  complex  structure  of  separation  ( S )  and  attachment  ( R )  lines.  The  initial  sepa¬ 
ration  lines  (S\  and  S2)  generated  by  the  individual  fin  shock  boundary  layer  interactions 
coalesce  to  form  a  node  at  the  centerline.  Corresponding  attachment  lines  (Ri  and  R2)  are 
located  near  the  fins.  Additional  separation  (S3  to  Sq)  and  attachment  lines  (R%  to  Rq)  form 
downstream.  The  central  point  ( S *)  is  a  saddle  point  in  the  experiment.  In  the  same  region, 
the  computation  of  Schmisseur  and  Gaitonde  (Fig.  93)  displays  a  node  (Afi),  while  the  com¬ 
putation  of  Panaras  (Fig.  94)  appears  in  closer  agreement  with  experiment.  Nonetheless, 
both  computed  surface  streamline  plots  show  overall  good  agreement  with  experiment. 

The  computed  and  experimental  surface  pressure  on  the  centerline  is  displayed  in  Fig.  95 
(Schmisseur  and  Gaitonde)  and  Fig.  96  (Panaras).  The  experimental  pressure  profile  displays 
a  rise  due  to  the  crossing  shock  interaction,  and  a  subsequent  drop  beginning  at  x  =  42 5  due 
to  the  expansion  originating  from  the  corners  C  in  Fig.  91.  The  subsequent  pressure  rise 
is  associated  with  the  reflection  of  the  shocks  from  the  fin  surfaces.  The  computed  surface 
pressure  of  Schmisseur  and  Gaitonde  shows  significant  deviation  from  experiment  beginning 
at  the  initial  plateau  in  surface  pressure.  The  computed  surface  pressure  of  Panaras  shows 
good  agreement  up  to  the  beginning  of  the  decrease  in  pressure  at  x/5  =  42.  However,  a 
significant  discrepancy  appears  thereafter. 

The  computed  and  experimental  surface  pressure  at  three  streamwise  locations  are  shown 
in  Figs.  97  to  102.  The  results  of  Schmisseur  and  Gaitonde  are  presented  in  Figs.  97,  99 
and  101,  and  the  results  of  Panaras  in  Figs.  98,  100  and  102.  Both  computations  show  good 
agreement  with  experiment  at  x  =  26.65  (Figs.  97  and  98).  This  location  is  upstream  of 
the  interaction  of  the  crossing  shocks,  and  thus  represents  individual  single  fin  interactions. 
At  x  =  37.15  (located  at  approximately  the  second  focal  point  N2  as  indicated  in  Fig.  93) 
both  computed  profiles  show  good  agreement  with  experiment,  except  for  a  overshoot  in 
centerline  pressure  in  Fig.  99  (Schmisseur  and  Gaitonde).  At  x  =  52.95  (see  Fig.  93),  both 
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computations  show  generally  good  agreement  with  experiment  except  for  an  overshoot  in 
pressure  at  the  centerline. 

In  summary,  the  following  conclusions  can  be  made: 

•  The  computed  surface  streamlines  are  in  good  agreement  with  experiment.  The  de¬ 
tailed  structure  of  separation,  attachment  and  foci  are  generally  accurately  reproduced 
in  the  simulations. 

•  The  computed  surface  pressure  shows  general  agreement  with  experiment  during  the 
initial  portion  of  the  interaction  and  off  centerline  throughout  the  entire  interaction. 
However,  significant  differences  between  computed  and  experimental  surface  pressure 
are  evident  on  the  centerline  of  the  interaction. 

5.4  Case  4:  Mach  5,  a\  =  23°,  «2  =  23° 

Schmisseur  and  Gaitonde  [200,  201,  205,  208]  and  Panaras  [227]  simulated  a  Mach  5,  = 

«2  =  23°  double  fin  at  Res  =  1.4  x  105  and  isothermal  wall  conditions  ( Tw/Taw  =  0.76). 
Experimental  data  of  Schiilein  and  Zheltovodov  [195,  207]  (an  additional  description  is  pre¬ 
sented  in  Zheltovodov  and  Maksimov  [228])  is  available  for  this  case.  The  flow  conditions 
are  shown  in  Table  20.  The  computations  were  performed  at  the  same  flow  conditions  as 
the  experiment.  The  experimental  configuration  is  shown  in  Fig.  103. 

Details  of  the  computations12  are  presented  in  Table  21.  Schmisseur  and  Gaitonde  utilized 
the  k—e  model  [220,  221]  with  low  Reynolds  number  and  compressibility  corrections.  Panaras 
[2]  employed  a  modified  Baldwin-Lomax  model. 

The  experimental  surface  oil  flow  is  presented  in  Fig.  104,  and  the  computed  surface 
streamlines  in  Figs.  105  (Schmisseur  and  Gaitonde)  and  106  (Panaras).  Good  agreement  is 
observed  between  the  computations  and  experiment. 

The  computed  and  experimental  surface  pressure  on  the  centerline  is  displayed  in  Fig.  107 
(Schmisseur  and  Gaitonde)  and  Fig.  108  (Panaras).  The  experimental  profile  is  similar  to 
the  18°  x  18°  configuration  except  for  the  absence  of  the  second  shock  reflection  on  the 
centerline  due  to  the  limited  length  of  the  flat  plate.  The  computed  profile  by  Schmisseur 
and  Gaitonde  shows  good  agreement  with  experiment  up  to  the  peak  experiment  pressure, 
but  significantly  overestimates  experiment  downstream.  The  computed  profile  by  Panaras 
also  shows  good  agreement  with  experiment  up  to  the  peak  pressure,  but  does  not  display 
the  expansion  region13  observed  in  the  experiment. 

12Schmisseur  and  Gaitonde  performed  a  grid  refinement  study  using  meshes  of  91  x  72  x  52,  181  x  143  x  103 
and  229  x  179  x  129  and  confirmed  that  the  finest  grid  solution  was  grid  converged. 

13A  separate  coarse  grid  computation  was  performed  by  Panaras  [227]  wherein  the  effect  of  the  finite  plate 
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The  computed  and  experimental  surface  pressure  at  three  streamwise  locations  are  shown 
in  Figs.  109  to  114.  The  results  of  Schmisseur  and  Gaitonde  are  shown  in  Figs.  109,  111 
and  113,  and  the  results  of  Panaras  in  Figs.  110,  112  and  114.  Both  computations  are  in 
close  agreement  with  the  available  experimental  data  at  x  —  22.15  which  is  upstream  of  the 
crossing  shock  interaction  (see  Figs.  107  and  108).  At  x  =  32.65,  located  in  the  midst  of 
the  crossing  shock  interaction,  both  computations  show  close  agreement  with  experiment, 
although  somewhat  overestimating  the  centerline  pressure.  At  x  =  48.45,  the  agreement  be¬ 
tween  computation  and  experiment  is  good  except  in  the  immediate  vicinity  of  the  centerline 
where  the  pressure  is  overestimated. 

In  summary,  the  following  conclusions  can  be  made: 

•  The  computed  surface  streamlines  display  good  agreement  with  experiment 

•  The  computed  surface  pressure  shows  general  agreement  with  experiment  during  the 
initial  portion  of  the  interaction  and  off  centerline  throughout  the  entire  interaction. 
However,  significant  differences  between  computed  and  experimental  surface  pressure 
are  evident  on  the  centerline,  similar  to  Case  3. 


6  Conclusions 

An  evaluation  of  the  CFD  capability  for  prediction  of  shock  wave  turbulent  boundary  layer 
interaction  was  performed  under  the  auspices  of  NATO  RTO  Working  Group  10.  Five 
separate  configurations  were  considered,  namely,  the  2-D  compression  corner,  2-D  shock 
impingement,  2-D  expansion-compression  corner,  3-D  single  fin  and  3-D  double  fin.  Recent 
DNS,  LES  and  RANS  simulations  were  solicited.  The  computed  results  were  compared  with 
experimental  data  where  available.  The  following  conclusions  may  be  drawn: 

•  DNS  and  LES 

—  DNS  and  LES  of  shock  wave  turbulent  boundary  layer  interaction  is  an  active 
area  of  research.  In  the  previous  review  of  shock  wave  turbulent  boundary  layer 
interaction  by  Knight  and  Degrez  [1]  published  in  1998,  there  were  no  DNS  or 
LES  results  available.  The  present  paper  includes  several  DNS  and  LES  results 
for  2-D  shock  wave  turbulent  boundary  layer  interactions. 

length  was  approximated  by  adding  a  divergent  section  to  the  flat  plate.  The  computed  centerline  pressure 
displayed  significantly  better  agreement  with  experiments,  but  the  spanwise  extent  of  the  central  structure 
was  smaller  than  in  the  experiment. 


RTO-TR-AVT-007-V3 


11  -49 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


ORGANIZATION 


—  A  comprehensive  assessment  of  the  capability  of  DNS  and  LES  was  not  possible 
for  two  reasons.  First,  the  DNS  and  LES  simulations  have  been  performed  at 
lower  Reynolds  numbers  Re, 5  than  the  experiment  for  all  cases  except  one  (be., 
the  2-D  shock  impingement).  Second,  the  DNS  and  LES  have  been  performed  for 
nominally  2-D  flows.  Both  reasons  are  attributable  to  the  computational  cost  of 
DNS  and  LES. 

—  Although  the  set  of  comparisons  between  DNS  and  LES  with  experiment  was  lim¬ 
ited,  nevertheless  it  is  clear  that  significant  progress  has  been  achieved  in  the  pre¬ 
diction  of  shock  wave  turbulent  boundary  layer  interaction  using  DNS  and  LES. 
In  several  cases,  the  computations  were  performed  at  Reynolds  numbers  close 
to  the  experiment  (and,  in  one  case,  equal  to  experiment),  and  good  agreement 
with  experimental  data  was  achieved.  Detailed  comparisons  included  surface  pres¬ 
sure,  skin  friction,  velocity,  temperature  and  turbulence  profiles.  However,  some 
significant  discrepancies  are  evident,  due  possibly  to  the  differences  in  Reynolds 
number. 

—  No  DNS  or  LES  results  for  surface  heat  transfer  in  shock  wave  turbulent  boundary 
layer  interaction  were  available.  This  is  particularly  disappointing,  since  prior 
RANS  simulations  of  strongly  separated  shock  wave  turbulent  boundary  layer 
interactions  have  failed  to  accurately  predict  heat  transfer  [1], 

•  RANS 

—  New  ideas  in  modeling  have  reinvigorated  RANS  simulations  for  shock  wave  tur¬ 
bulent  boundary  layer  interaction.  These  new  ideas  include  realizability  (weak 
nonlinearity)  and  specific  physical  models  which  incorporate  known  flowfield  be¬ 
havior. 

—  The  weakly  nonlinear  Wilcox-Durbin  model  successfully  predicts  secondary  sep¬ 
aration  for  the  3-D  single  bn  in  agreement  with  experiment. 

—  The  linear  and  weakly  nonlinear  Wilcox-based  models  fail  to  accurately  predict 
surface  heat  transfer  for  the  separated  3-D  double  bn. 

—  The  surface  streamline  pattern  for  the  3-D  double  bn  Case  Nos.  3  and  4  is  pre¬ 
dicted  by  both  the  conventional  k  —  e  and  Baldwin- Lomax-Panaras  models. 

—  The  surface  pressure  for  the  3-D  double  bn  is  accurately  predicted  by  the  conven¬ 
tional  k  —  e  and  Baldwin-Lomax-Panaras  models  only  in  the  initial  region  of  the 
interaction.  Signihcant  diberences  between  the  computational  and  experimental 
centerline  pressure  are  evident  downstream. 
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7  Future  Needs 

The  following  research  needs  for  shock  wave  turbulent  boundary  layer  interaction  may  be 
identified: 

•  DNS  and  LES  simulations  are  needed  at  Reynolds  numbers  corresponding  to  the  ex¬ 
periment.  This  is  essential  to  end  speculation  regarding  the  reasons  for  discrepancy 
between  simulation  and  experiment. 

•  DNS  and  LES  simulations  are  needed  for  3-D  configurations.  All  DNS  and  LES  results 
to  date  have  been  for  nominally  2-D  flows. 

•  DNS  and  LES  simulations  are  needed  to  compare  with  experimental  data  for  heat  trans¬ 
fer.  Neither  DNS,  LES  nor  RANS  models  to  date  have  successfully  predicted  surface 
heat  transfer  in  strongly  separated  shock  wave  turbulent  boundary  layer  interactions. 

•  RANS  modeling  needs  continued  emphasis  on  weakly  nonlinear  corrections  to  two 
equation  models,  and  continued  development  of  configuration-specific  modifications  to 
simple  models. 
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Figure  1:  2-D  compression  corner 
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Figure  3:  2-D  expansion-compression  corner 
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Figure  4:  3-D  single  fin 
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Figure  5:  3-D  double  fin 
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Figure  6:  Surface  pressure  for  2-D  compression  corner  and  shock  impingement  at  =  2.96 
and  a  =  25°  [19] 
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Figure  7:  Separation  length  vs  Res  [24] 
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Figure  8:  Turbulence  structure  in  2D  compression  corner 
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Figure  9:  Separation  length  for  RANS 
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Figure  10:  Instantaneous  velocity  and  shock 
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Figure  11:  Mean  wall  pressure 
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Figure  12:  Mean  velocity 
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Figure  13:  Mean  static  temperature 
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Figure  14:  Reynolds  streamwise  stress 
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Figure  15:  Reynolds  shear  stress 
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Figure  16:  Generation  of  inflow  conditions 
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Figure  17:  Inflow  mean  velocity 
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Figure  18:  Inflow  rms  streamwise  velocity 
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Figure  19:  Mean  surface  pressure 
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Figure  20:  Mean  skin  friction 
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Figure  21:  Mean  velocity 
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Figure  22:  Mean  static  temperature 


RTO-TR-AVT-007-V3 


11  -93 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


ORGANIZATION 


DNS 


Figure  23:  Spanwise  TKE  spectra  at  y+  =  7.3 (y/5  =  0.034) 


11  -94 


RTO-TR-AVT-007-V3 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


DNS 


Figure  24:  Spanwise  TKE  spectra  at  y+  =  132 (y/5  =  0.61) 
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Figure  25:  Turbulence  statistics  (-  DNS,. . .  LES  (ADM)) 
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Figure  26:  Separation  length  for  DNS  and  LES 
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Figure  27:  Inflow  velocity 
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Figure  28:  Inflow  Reynolds  streamwise  stress 
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Figure  29:  Inflow  Reynolds  shear  stress 
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Figure  30:  Mean  wall  pressure 
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Figure  31:  Standard  deviation  of  wall  pressure 
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Figure  32:  Skewness  of  wall  pressure 
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Figure  33:  Time  history  of  wall  pressure 
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Figure  34:  Mean  skin  friction 
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Figure  35:  Velocity  profiles 
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Figure  36:  Mean  wall  pressure 
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Figure  37:  Mean  skin  friction 
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Figure  38:  Mean  wall  pressure 
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Figure  39:  Mean  skin  friction  coefficient 
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Figure  41:  Instantaneous  Schlieren  image 
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Figure  42:  Velocity  (A  — ,  B  A,  C - ,  D  o,  HWA  o) 
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Figure  43:  Skin  friction  (A  — ,  B  A,  C - ,  Do,  HWA  o) 
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Figure  44:  Displacement  thickness  (A  — ,  B  A,  C - ,  D  o,  HWA  o) 


RTO-TR-AVT-007-V3 


11  - 115 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


ORGANIZATION 


Figure  45:  Velocity  (A  (x  =  351  mm)  — ,  A  (x  —  345  mm)  —  •  — ,  B  A,  C - ,  D  o,  HWA 

o,  LDA  □) 
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Figure  46:  Rms  velocity  (A  (x  =  326  mm)  — ,  A  (x  —  320  mm)  —  •  — ,  B  A,  C - ,  D  o, 

LDA  □) 
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Figure  47:  Rms  velocity  (A  (x  =  351  mm)  — ,  A  (x  =  345  mm)  —  •  — ,  B  A,  C - ,  D  o, 

HWA  o,  LDA  □) 
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x  =  125  mm 


x  =  0  mm 


Figure  48:  Instantaneous  Schlieren  images  of  (a)  the  flat-plate  boundary  layer,  (b)  7°  centered 
expansion,  and  (c)  14°  centered  expansion 
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[b) 


Figure  49:  Instantaneous  FRS  visualization  of  (a)  the  7°  centered  and  (b)  gradual  expansion 
regions.  The  white  lines  indicate  the  position  of  the  model  surface. 
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Figure  50:  Separation  length 
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Figure  51:  Mean  static  pressure 
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Figure  52:  Mean  streamlines 
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Figure  53:  Mean  velocity 
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Figure  54:  Surface  pressure 
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Figure  55:  Skin  friction  coefficient 
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Figure  56:  Flow  model  proposed  by  Kubota  and  Stollery[146] 
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Figure  57:  Flow  model  of  Alvi  and  Settles[148] 
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Cross-section  of  the  vortices 


Figure  58:  Perspective  view  of  the  conical  separation  vortices  and  of  the  shock  waves  (Panaras 
[150] 
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Figure  59:  Conical  projections  of  sections  of  the  flow  on  the  outflow  plane  (Panaras[150]) 
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Figure  60:  Skin  friction  lines  and  cross-section  of  the  vortices  (Panaras  and  Stanewsky[155]) 
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Figure  61:  3-D  single  fin  regimes  (Zheltovodov  et  a/[156]) 
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(a)  Cross-section  of  the  flow  where  the  thick  line  denotes  the  y  =  0.78  stream  line 


Vortex  Core 


Development  of  the  y=0  ' 
stream  surface 


(b)  Perspective  development  of  the  y  =  0.75  stream  surface 


Figure  62:  Structure  of  the  separated  flow  in  a  swept  shock/turbulent  boundary-layer  inter¬ 
action 


RTO-TR-AVT-007-V3 


11  - 133 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


ORGANIZATION 


4.0 

3.0 

5$°  2.0 

>■, 

1.0 
0.0 

Z/5o 

(a)  Moo  =  4.0,  a  =  20° 

4.0 

3.0 

S$°  2.0 

>. 

1.0 
0.0 

0  2  4  6  8  10  12  14  16  18  20  22  24 

z/s0 

(b)  Moo  =  3.0,  a  =  10° 

Figure  63:  Cross-sections  of  various  flows  (Panaras  [2]) 
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Figure  64:  Grid  for  3-D  single  fin 
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Figure  66:  Surface  streamlines  for  Case  1  (exper) 
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Figure  67:  Surface  streamlines  for  Case  1  (WD+) 
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Figure  68:  Surface  streamlines  for  Case  3  (WI) 
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Figure  69:  Surface  streamlines  for  Case  3  (exper) 


11  -140 


RTO-TR-AVT-007-V3 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  interactions 


Z1  (mm) 


Figure  70:  Surface  streamlines  for  Case  3  (WD+) 
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Figure  71:  Surface  pressure  for  Case  1 
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P/P,  (R  =  88.9  mm) 


Figure  72:  Surface  pressure  for  Case  2 
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Figure  73:  Surface  pressure  for  Case  3 


11  -144 


RTO-TR-AVT-007-V3 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


Cfx1000  (R  =  88.9  mm) 


Figure  74:  Skin  friction  for  Case  2 
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Figure  75:  Turbulence  kinetic  energy  for  Case  3  (WI) 
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Figure  76:  Turbulence  kinetic  energy  for  Case  3  (WD+) 
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Figure  77:  3-D  Double  Fin  Flow  Structure  [177,  173,  178] 
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Figure  78:  3-D  Double  Fin  Flow  Structure  [189,  192,  193,  197] 
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Figure  79:  3-D  Double  Fin  Flow  Structure  [189,  192,  193,  197] 
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Figure  80:  3-D  Double  Fin  Flow  Structure  [189,  192,  193,  197] 
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Figure  81:  3-D  Double  Fin  Flow  Structure  [159] 
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c)  Numerical  solutions  in  the  central  separation  region 


Figure  82:  Topology,  skin  friction  and  heat  transfer 
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Figure  83:  Surface  pressure  and  heat  transfer  (an  =  15°,  a?2  =  15°) 
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Figure  84:  3-D  double  fin  flow  structure 


RTO-TR-AVT-007-V3 


11  - 155 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


ORGANIZATION 


Ch  x  1 03  :  0.67  1.00  2.00  2.70  3.00  4.00 


Figure  85:  Computed  surface  heat  transfer 
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Figure  86:  Surface  flow  visualization 
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Figure  87:  Surface  pressure  and  heat  transfer  (an  =  7°,a2  =  11°) 
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Figure  88:  Surface  pressure  (ctu  =  7°,  «2  =  11°) 
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Figure  89:  Surface  pressure  (ctu  =  7°,  «2  =  11°) 


11 -160 


RTO-TR-AVT-007-V3 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


N  =  0  1  2  3  4  5 

zr=  2.6  2.2  1.8  0.75  0.33  0.07 


'/l,og  (+N) 


Figure  90:  Turbulence  length  scale  ratio  (au  =  70,cc2  =  11°) 
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Figure  91:  Double  fin  configuration 
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Figure  92:  Experimental  surface  oil  visualization 
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spanwise  pressure  survey  locations 


Figure  93:  Surface  streamlines 
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Figure  94:  Surface  streamlines 
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Figure  95:  Surface  pressure  on  centerline 
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Figure  96:  Surface  pressure  on  centerline 
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Figure  97:  Surface  pressure  at  x  —  26.65 
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Figure  98:  Surface  pressure  at  x  —  26.65 


RTO-TR-AVT-007-V3 


11  - 169 


Advances  in  CFD  Prediction  of  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 


ORGANIZATION 


-10  -5  0 


Z(8) 

Figure  99:  Surface  pressure  at  x  —  37.15 
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Figure  100:  Surface  pressure  at  x  —  37.15 
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Figure  101:  Surface  pressure  at  x  —  52.96 
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Figure  102:  Surface  pressure  at  x  —  52.95 
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Figure  103:  Double  fin  configuration 
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Figure  104:  Experimental  surface  oil  visualization 
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Figure  105:  Surface  streamlines 
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Figure  106:  Surface  streamlines 
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Figure  107:  Surface  pressure  on  centerline 
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Figure  108:  Surface  pressure  on  centerline 
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Figure  109:  Surface  pressure  at  x  —  22.15 
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Figure  110:  Surface  pressure  at  x  —  22.1  S 
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Figure  111:  Surface  pressure  at  x  —  32.65 
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Figure  112:  Surface  pressure  at  x  —  32.6 5 
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Figure  113:  Surface  pressure  at  x  —  48.4 5 
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Figure  114:  Surface  pressure  at  x  —  48.45 
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Table  1:  RANS  calculations  of  2-D  Compression  Corner 


Number 

Ref. 

Model 

a 

Mach 

Res  x  10"6 

Lsep/ ^ 

Lsep/ Lc 

1 

Shang  [60] 

Equilibrium 

25° 

2.96 

0.1 

5.0 

15.3 

2 

Shang  [60] 

Frozen 

25° 

2.96 

0.1 

13.0 

39.7 

3 

Shang  [60] 

Relaxation 

25° 

2.96 

0.1 

10 

30.5 

4 

Horstman  [61] 

Baseline 

20° 

2.85 

1.65 

1.45 

8.7 

5 

Horstman  [61] 

Relaxation 

20° 

2.85 

1.65 

2.69 

16.1 

6 

Horstman  [61] 

Pressure  Gradient 

20° 

2.85 

1.65 

0.81 

4.87 

7 

Horstman  [61] 

Relaxation+ 
Pressure  Gradient 

20° 

2.85 

1.65 

1.62 

9.75 

8 

Horstman  [61] 

Baseline 

24° 

2.85 

1.33 

4.0 

14.5 

9 

Horstman  [61] 

Relaxation 

24° 

2.85 

1.33 

4.7 

17.1 

10 

Horstman  [61] 

Pressure  Gradient 

24° 

2.85 

1.33 

2.74 

9.96 

11 

Horstman  [61] 

Relaxation+ 
Pressure  Gradient 

24° 

2.85 

1.33 

3.6 

13.1 

12 

Visbal, Knight  [62] 

B-L 

16° 

2.9 

1.6 

1.0 

14.4 

13 

Visbal, Knight  [62] 

Modihed  B-L 

16° 

2.9 

1.6 

0.2 

2.89 

14 

Visbal, Knight  [62] 

Relaxation 

16° 

2.9 

1.6 

0.6 

8.66 

15 

Visbal, Knight  [62] 

Constant  Yrnax 

16° 

2.9 

1.6 

0.2 

2.89 

16 

Visbal, Knight  [62] 

B-L 

20° 

2.9 

1.6 

1.6 

11.5 

17 

Visbal, Knight  [62] 

Modihed  B-L 

20° 

2.9 

1.6 

0.9 

6.46 

18 

Visbal, Knight  [62] 

Relaxation 

20° 

2.9 

1.6 

1.8 

12.9 

19 

Visbal, Knight  [62] 

Modihed  B-L 

24° 

2.9 

1.6 

2.05 

7.25 

20 

Visbal, Knight  [62] 

Relaxation 

24° 

2.9 

1.6 

2.95 

10.4 

21 

Ong, Knight  [63] 

B-L 

16° 

1.96 

0.25 

2.15 

10.1 

22 

Ong, Knight  [63] 

B-L 

16° 

2.83 

1.6 

0.6 

2.81 

23 

Ong, Knight  [63] 

B-L 

20° 

2.83 

1.6 

1.75 

8.19 

24 

Horstman  [64] 

k  —  e 

25° 

3.0 

0.133 

4.89 

16.3 

25 

Borisov  et  al  [65] 

k  —  uj 

25° 

3.0 

0.139 

7.38 

24.8 
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Table  2:  Flow  Conditions 


Reference 

Data 

Moo 

Res  x  10“4 

Urbin  et  al  [85] 

LES 

3.0 

2.0 

Kannepalli  et  al  [86] 

LES 

2.88 

2.0 

Zheltovodov  et  al[ 35,  36,  37] 

E 

2.9 

7.5  -  11 
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Table  3:  Details  of  Computations 


Reference  Type 

Urbin  et  al  [85]  LES 

Kannepalli  et  al  [86]  LES 


Cells 

SA 

TA 

1.6 

x  106 

2nd 

2nd 

300  x 

65  x  56 

4th 

2nd/4th 
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Table  4:  Flow  Conditions 


Reference 

Data 

Moo 

Res  x  10 

Adams  [90,  91,  44] 

DNS 

3 

2.1 

Rizzetta  et  al  [92] 

DNS,  LES 

3 

2.1 

Stolz  et  al  [93] 

LES 

3 

2.1 
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Table  5:  Details  of  Computations 


Reference 

Type 

Adams  [90,  91,  44] 

DNS 

Rizzetta  et  al{ 92] 

DNS 

LES-S 

LES-D 

Stolz  et  al  [93] 

LES 

Grid 

SA 

TA 

1000x180x80 

5th 

3rd 

421x151x81 

6th 

2nd 

421x151x81 

6th 

2nd 

421x151x81 

6th 

2nd 

334  x  91  x  31 

6th 

3rd 
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Table  6:  Flow  Conditions 


Reference 

Data 

Res  x  10~6 

Rizzetta  et  al[ 97] 

LES 

3.0 

.021 

Kannepalli  et  al[ 86] 

LES 

2.88 

0.02 

Settles  et  al  [98] 

E 

2.84-2.87 

1.5-  1.6 

Dolling  et  al  [49] 

E 

2.90-2.95 

0.81  -  1.4 

Smits  et  al  [99] 

E 

2.79-2.87 

1.6 

Zheltovodov  et  al[ 3,  35,  36,  37,  100] 

E 

2.88 

0.063 
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Table  7:  Details  of  Computations 


Reference 

Type 

Grid 

SA 

TA 

Rizzetta  et  al[ 97] 

LES 

421  x  151  x  181 

6th 

2nd 

Kannepalli  et  al  [86] 

LES 

300  x  65  x  56 

4th 

2nd/4th 
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Table  8:  Flow  Conditions 
Reference  Data 

Yan  et  al[ 59]  LES 

Zheltovodov  et  al[ 3,  35,  36,  37,  100]  E 


Moo 

2.88 

2.88 


Res  x  10“4 
2.0 
6.4 
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Table  9: 
Reference 
Yan  et  al[ 59] 


ORGANIZATION 


Details  of  Computations 
Type  Cells  SA  TA 

LES  2.0  x  106  2nd  2nd 
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Table  10:  Flow  Conditions 
Reference  Data 

Gamier  et  aZ[101,  102]  LES  2.3 
Dcleuze[103],  Laurent[104]  E  2.3 


Res  x  10  4 
6.0 
6.0 
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Table  11:  Details  of  Computations 


Reference 

Type 

Garnier[101,  102] 

LES-X 

(a) 

LES-X 

(b) 

LES-X 

(c) 

LES-M 

(d) 

Grid 

SA 

TA 

255  x  55  x  151 

2nd 

2nd 

255  x  110  x  151 

2nd 

2nd 

510  x  55  x  151 

2nd 

2nd 

510  x  55  x  151 

2nd 

2nd 
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Table  12:  Flow  Conditions 
Reference  Data 

Knight  et  a/[139]  LES 

Zheltovodov[23,  24,  37,  39,  40,  138,  136,  41]  E 


Moo 

2.88 

2.88 


Res  x  10  4 

2.0 

4.4-50. 
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Table  13:  Details  of 
Reference  Type 
Knight  et  a/[139]  LES 


Computations 
Cells  SA 
2.4  x  106  2nd 


TA 

2nd 
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Table  14:  Flow  Conditions 


Case 

Reference 

Data 

a 

Res 

xl0~ 

1 

Thivet[161] 

RANS 

3.0 

15° 

1.9 

Zheltovodov[154,  143] 

E 

3.0 

15° 

1.9 

2 

Thivet[161] 

RANS 

4.0 

20° 

2.1 

Kim  et  al  [158] 

E 

4.0 

20° 

2.1 

3 

Thivet[161] 

RANS 

4.0 

30.6° 

1.6 

Zheltovodov  et  al  [154] 

E 

4.0 

30.6° 

1.6 
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Table  15:  Details  of  Computations 


Case 

Reference 

Type 

Grid 

SA 

TA 

1 

Thivet[161] 

RANS 

128  x  160  x  160 

2nd 

SS 

2 

Thivet[161] 

RANS 

128  x  160  x  160 

2nd 

ss 

3 

Thivet[161] 

RANS 

112  x  160  x  160 

2nd 

SS 
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Table  16:  Flow  Conditions 


Case 

Reference 

Data 

Moo 

«2 

Res 

xl0“ 

1 

Thivet[164,  202,  204,  209] 

RANS 

3.92 

15° 

15° 

3.0 

Zheltovodov  et  aZ[187,  194,  196, 

216] 

E 

3.92 

15° 

15° 

3.0 

2 

Thivet[164,  202,  204,  209] 

RANS 

3.92 

7° 

11° 

3.0 

Zheltovodov  et  a/[187,  194,  196, 

216] 

E 

3.92 

7° 

11° 

3.0 
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Table  17:  Details  of  Computations 


Case 

Reference 

Type 

Grid 

SA 

TA 

1 

Thivet[164] 

RANS 

192  x  200  x  104 

2nd 

SS 

2 

Thivet[164] 

RANS 

96  x  80  x  160 

2nd 

ss 
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Table  18: 

Reference 

Schmisseur  et  a/[200] 

Panaras[227] 

Schiilein  and  Zheltovodov[195,  228,  207] 


Conditions 


Data 

Moo 

Ot  i 

OL  2 

Res 

RANS 

4.96 

o 

00 
T— 1 

o 

00 
T— 1 

1.39  x  105 

RANS 

4.96 

00 

o 

o 

00 
T— 1 

1.39  x  105 

E 

4.96 

00 

o 

o 

00 
T— 1 

1.39  x  105 

Flow 
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Table  19: 

Reference 

Schmisseur  et  a/[200] 
Panaras[227] 


Details  of  Computations 


Type 

Grid 

SA 

TA 

RANS 

153  x  143  x  103 

2nd 

SS 

RANS 

127  x  101  x  91 

2nd 

ss 
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Table  20: 

Reference 

Schmisseur  et  a/[200] 

Panaras[227] 

Schiilein  and  Zheltovodov[195,  228,  207] 


Conditions 


Data 

Moo 

c%i 

Qi2 

Res 

RANS 

4.96 

23° 

23° 

1.39  x  105 

RANS 

4.96 

23° 

23° 

1.39  x  105 

E 

4.96 

23° 

23° 

1.39  x  105 

Flow 
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Table  21: 

Reference 

Schmisseur  et  a/[200] 
Panaras[227] 


Details  of  Computations 


Type 

Grid 

SA 

TA 

RANS 

229  x  179  x  129 

2nd 

SS 

RANS 

117  x  101  x  91 

2nd 

ss 
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